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We review and analyze the available information for nuclear fusion cross 



sections that are most important for solar energy generation and solar neu- 
trino production. We provide best values for the low-energy cross-section 
factors and, wherever possible, estimates of the uncertainties. We also de- 
scribe the most important experiments and calculations that are required in 
order to improve our knowledge of solar fusion rates. 
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I. INTRODUCTION 



This section describes in Sec. !I A] the reasons why a critical analysis of what is known 
about solar fusion reactions is timely and important, summarizes in Sec. fB| the process by 



which this collective manuscript was written , and provides in Sec. !1 C! a brief outline of the 
structure of the paper. 
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A. Motivation 



The original motivation of solar neutrino experiments was to use the neutrinos "..to 
see into the interior of a star and thus verify directly the hypothesis of nuclear energy 
generation in stars" (Bahcall, 1964; Davis, 1964). This goal has now been achieved by 
four pioneering experiments: Homestake (Davis, 1994), Kamiokande (Fukuda et al, 1996), 
GALLEX (Kirsten et al, 1997), and SAGE (Gavrin et al, 1997). These experiments provide 
direct evidence that the stars shine and evolve as the result of nuclear fusion reactions among 
light elements in their interiors. 

Stimulated in large part by the precision obtainable in solar neutrino experiments and 
by solar neutrino calculations with standard models of the sun, our knowledge of the low- 
energy cross sections for fusion reactions among light elements has been greatly refined by 
many hundreds of careful studies of the rates of these reactions. The rate of progress was 
particularly dramatic in the first few years following the proposal of the chlorine (Homestake) 
experiment in 1964. 

In 1964, when the chlorine solar neutrino experiment was proposed (Davis, 1964; Bahcall, 
1964), the rate of the ^He-^He reaction was estimated (Good, Kunz, and Moak, 1954; Parker, 
Bahcall, and Fowler, 1964) to be 5 times slower than the current best estimate and the 
uncertainty in the low-energy cross section was estimated (Parker, Bahcall, and Fowler, 
1964) to be "as much as a factor of 5 or 10." Since the ^He-^He reaction competes with 
the ^He-^He reaction — which leads to high energy neutrinos — the calculated fluxes for the 
higher energy neutrinos were overestimated in the earliest days of solar neutrino research. 
The most significant uncertainties, in the rates of the ^He-^He, the ^He-^He, and the ^Be-p 
reactions, were much reduced after just a few years of intensive experimental research in the 
middle and late 1960s (Bahcall and Davis, 1982). 

Over the past three decades, steady and impressive progress has been made in refining 
the rates of these and other reactions that produce solar energy and solar neutrinos. (For 
reviews of previous work on this subject, see, e.g.. Fowler, Caughlan, and Zimmerman, 1967, 
1975; Bahcall and Davis, 1982; Clayton, 1983; Fowler, 1984; Parker, 1986; Rolfs and Rodney, 
1988; Caughlan and Fowler, 1988; Bahcall and Pinsonneault, 1992, 1995; Parker and Rolfs, 
1991). An independent assessment of nuclear fusion reaction rates is being conducted 
by the European Nuclear Astrophysics Compilation of Reaction Rates (NACRE) (see, e.g., 
Angulo, 1997); the results from this compilation, which has broader goals than our study 
and in particular does not focus on precision solar rates, are not yet available. 

However, an unexpected development has occurred. The accuracy of the solar neutrino 
experiments and the precision of the theoretical predictions based upon standard solar mod- 
els and standard electroweak theory have made possible extraordinarily sensitive tests of new 
physics, of physics beyond the minimal standard electroweak model. Even more surprising 
is the fact that, for the past three decades, the neutrino experiments have consistently dis- 
agreed with standard predictions, despite concerted efforts by many physicists, chemists, 
astronomers, and engineers to find ways out of this dilemma. 

The four pioneering solar neutrino experiments together provide evidence for physics 
beyond the standard electroweak theory. The Kamiokande (Fukuda et al, 1996) and the 
chlorine (Davis, 1994) experiments appear to be inconsistent with each other if nothing 
happens to the neutrinos after they are created in the center of the sun (Bahcall and Bethe, 
1990). Moreover, the well calibrated gallium solar neutrino experiments GALLEX (Kirsten 
et al, 1997), and SAGE (Gavrin et al, 1997) are interpreted, if neutrinos do not oscillate 
or otherwise change their states on the way to the earth from the solar core, as indicating 
an almost complete absence of ^Be neutrinos. However, we know [see discussion of Eq. (pSj) 
in Sec. [VII|| that the ''Be neutrinos must be present, if there is no new electroweak physics 
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occurring, because of the demonstration that neutrinos are observed by the Kamiokande 
solar neutrino experiment. Both ''Be and ^B neutrinos are produced by capture on ^Be ions. 

New solar neutrino experiments are currently underway to test for evidence of new physics 
with exquisitely precise and sensitive techniques. These experiments include a huge pure 
water Cerenkov detector known as Super-Kamiokande (Suzuki, 1994; Totsuka, 1996), a 
kiloton of heavy water, SNO, that will study both neutral and charged currents (Ewan et 
ai, 1987, 1989; McDonald, 1995), a large organic scintillator, BOREXINO, that will in- 
vestigate lower energy neutrinos than has previously been possible (Arpesella et al., 1992; 
Raghavan, 1995), and a 600 ton liquid argon time projection chamber, ICARUS, that will 
provide detailed information on the surviving '^B flux (Rubbia, 1996; ICARUS collab- 
oration, 1995; Bahcall et ai, 1986). With these new detectors, it will be possible to 
search for evidence of new physics that is independent of details of solar model predic- 
tions. [Discussions of solar neutrino experiments and the related physics and astronomy 
can be found at, for example, ^ttp: //www. hep. anl.gov/NDK/Hypertext /nuindustry.html] , 



http: / /neutrino. pc.helsinki.fi/neutrino/ , and http:/ /www.sns.ias.edu^jnbl . 



However, our ability to interpret the existing and new solar neutrino experiments is 
limited by the imprecision in our knowledge of the relevant nuclear fusion cross sections. To 
cite the most important example, the calculated rate of events in the Super-Kamiokande 
and SNO solar neutrino experiments is directly proportional to the rate measured in the 
laboratory at low energies for the '^Be(p, 7)^B reaction. This reaction is so rare in the sun, 
that the assumed rate of ^Be(p, 7)^B has only a negligible effect on solar models and therefore 
on the structure of the sun. The predicted rate of neutrino events in the interval 2 MeV to 
15 MeV is directly proportional to the measured laboratory rate of the '^Be(p, 7)^B reaction. 
Unfortunately, the low-energy cross-section factor for the production of ^B is the least well 
known of the important cross sections in the pp chain. 

We will concentrate in this review on the low-energy cross section factors, S, that deter- 
mine the rates for the most important solar fusion reactions. The local rate of a non-resonant 
fusion reaction can be written in the following form (see, e.g., Bahcall, 1989): 



{av) = 1.3005 X 10 



-15 



Z1Z2 

'at! 



1/3 

/S'eflf exp (-r) cm^ s"^ (1) 



Here Zi, Z2 are the nuclear charges of the fusing ions, Ai,A2 are the atomic mass numbers, 
A the reduced mass AiA2/{Ai + A2), Tg is the temperature in units of 10^ K, and the cross- 
section factor S'eff (defined below) is in keV b. The most probable energy, Eq, at which the 
reaction occurs is 



En 



{'naZiZ2kTf{mAc^ /2)^^^ 
1.22QA{ZlZlAT^f'^ keV. (2) 



The energy Eq is also known as the Gamow energy. The exponent r that occurs in Eq. (|I]) 
dominates the temperature dependence of the reaction rate and is given by 

r = 3Eo/A;T = 42.487 (zlZlAT^^) ^'^ . (3) 



For all the important reactions of interest in solar fusion, r is in the range 15 to 40. The 
quantity / is a correction factor due to screening first calculated by Salpeter (1954) and 
discussed in this paper in Sec. |11 C] . The quantity Seff is the effective cross section factor 
for the fusion reaction of interest and is evaluated at the most probable interaction energy, 
Eq. To first order in (Bahcall, 1966), 
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S,s = S{Eo) 1 + r- 



5 5S'Eo S"El 

\ H 

12 25 S 



(4) 



E=Eo. 



Here S' = dS/dE. In most analyses in the literature, the values of S and associated 
derivatives are quoted at zero energy, not at Eq. In order to relate to the usual formulae, 
one must express the relevant quantities in terms of their values a.t E = 0. The appropriate 
connection is 



'S'efr(-E'o) 
5(0) [1 



5 S' {Eo + %kT) s"Eo 
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S 



Eo 89 , ^ 
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(5) 



In some contexts, Sc{{{Eq) is referred to as simply the 'S'-factor' or 'the low-energy ^-factor'. 

For standard solar models (cf. Bahcall, 1989), the fusion energy and the pp neutrino 
flux are generated over a rather wide range of temperatures, 8 < < 16. The other 
important fusion reactions and neutrino fluxes are generated over a more narrow range of 
physical conditions. The neutrino flux is created in the most restricted temperature 
range, 13 < Tg < 16 . The mass density (in g cm~^) is given approximately by the relation 
p = 0.04T| in the temperature range of interest. 

The approximate dependences of the solar neutrino fluxes on the different low-energy 
nuclear cross-section factors can be calculated for standard solar models. The most impor- 
tant fluxes for solar neutrino experiments that have been carried out so far, or which are 
currently being constructed, are the low energy neutrinos from the fundamental pp reac- 
tion, 0(pp), the intermediate energy '^Be line neutrinos, 0(^Be), and the rare high-energy 
neutrinos from ^B decay, 0(^B). The pp neutrinos are the most abundant experimentally- 
accessible solar neutrinos and the ^B neutrinos have the smallest detectable flux, according 
to the predictions of standard models (Bahcall, 1989). 

Let ^ii, 533, and 5*34 be the low-energy, nuclear cross-section factors (defined in Sec. P A|) 
for the pp, ^He + ^He, and ^He + ^He reactions and let 5*17 and Se- 7 be the cross-section 
factors for the capture by ^Be of, respectively, protons and electrons. Then (Bahcall 1989) 



^34 



(6a) 



(/.(^Be)oc5rr''53-3°-^''^°-'' 



■^34 



and 



(6b) 



^B)cx5ri^-«53-3°-^° 



c0.81 cl.O c-1.0 



(6c) 



Nuclear fusion reactions among light elements both generate solar energy and produce 
solar neutrinos. Therefore, the observed solar luminosity places a strong constraint on 
the current rate of solar neutrino generation calculated with standard solar models. In 
addition, the shape of the neutrino energy spectrum from each neutrino source is unaffected, 
to experimental accuracy, by the solar environment. A good fit to the results from current 
solar neutrino experiments is not possible, independent of other, more model-dependent solar 
issues provided nothing happens to the neutrinos after they are created in the sun (see, e.g., 
Castellani et al. 1997; Heeger and Robertson, 1996; Bahcall, 1996; Hata, Bludman, and 
Langacker, 1994, and references therein). 
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But, the ultimate limit of our ability to extract astronomical information and to infer 
neutrino parameters will be constrained by our knowledge of the spectrum of neutrinos 
created in the center of the sun. Returning to the example of the neutrinos, the total 
flux (independent of flavor) of these neutrinos will be measured in the neutral current ex- 
periment of SNO, and-using the charged current measurements of SNO and ICARUS-in 
Super-Kamiokande. This total flux is very sensitive to temperature, 0(^B) ~ S^T^'^ (Bah- 
call and Ulmer, 1996), where T is the central temperature of the sun. Therefore, our abihty 
to test solar model calculations of the central temperature profile of the sun is limited by 
our knowledge of 5*17. 

Existing or planned solar neutrino experiments are expected to determine whether the 
energy spectrum of electron type neutrinos created in the center of the sun is modified 
by physics beyond standard electroweak theory. Moreover, these experiments have the 
capability of determining the mechanism, if any, by which new physics is manifested in 
solar neutrino experiments and thereby determining how the original neutrino spectrum is 
altered by the new physics. Once we reach this stage, solar neutrino experiments will provide 
precision tests of solar model predictions for the rates at which nuclear reactions occur in 
the sun. 

After the neutrino physics is understood, neutrino experiments will determine the average 
ratio in the solar interior of the ^He-^He reaction rate to the rate of the '^He-^He reaction. 
This solar ratio of reaction rates, -R33/-R34, can be inferred directly from the measured total 
fiux of ^Be and pp neutrinos (Bahcall, 1989). The comparison of the measured and the 
calculated ratio of -R33/-R34 will constitute a stringent and informative test of the theory of 
stellar interiors and nuclear energy generation. In order to extract the inherent information 
about the solar interior from the measured ratio, we must know the nuclear fusion cross 
sections that determine the branching ratios among the different reactions in the pp chain. 



B. The origin of this work 

This paper originated from our joint efforts to critically assess the state of the nuclear 
physics important to the solar neutrino problem. There are two motivations for taking on 
such a task at this time. First, we have entered a period where the sun, and solar mod- 
els, can be probed with unprecedented precision through neutrino flux measurements and 
helioseismology. It is therefore important to assess how uncertainties in our understanding 
of the underlying nuclear physics might affect our interpretation of such precise measure- 
ments. Second, as the importance of the solar neutrino problem to particle physics and 
astrophysics has grown, so has also the size of the community interested in this problem. 
Many of the interested physicists are unfamiliar with the decades of effort that have been 
invested in extracting the needed nuclear reaction cross sections, and thus uncertain about 
the quality of the results. The second goal of this paper is to provide a critical assessment 
of the current state of solar fusion research, describing what is known while also delineating 
the possibilities for further reducing nuclear cross-section uncertainties. 

In order to achieve these goals, an international collection of experts on nuclear physics 
and solar fusion — representing every speciality (experimental and theoretical) and every 
point of view (often conflicting points of view) — met in a workshop on "Solar Fusion Reac- 
tions." In particular, the participants included experts on all the major controversial issues 
discussed in widely circulated preprints or in the published literature. The workshop was 
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held at the Institute for Nuclear Theory, University of Washington, February 17-20, 1997.Q 
The goal of the workshop was to initiate critical discussions evaluating all of the existing 
measurements and calculations relating to solar fusion and to recommend a set of standard 
parameters and their associated uncertainties on which all of the participants could agree. 
To achieve this goal, we undertook ab initio analyses of each of the important solar fusion 
reactions; previously cited reviews largely concentrated on incremental improvements on 
earlier work. This paper is our joint work and represents the planned culmination of the 
workshop activities. 

At the workshop, we held plenary sessions on each of the important reactions and also 
intensive specialized discussions in smaller groups. The discussions were led by the fol- 
lowing individuals: extrapolations (K. Langanke), electron screening (S. Koonin), pp (M. 
Kamionkowski), ^He +^ He (C. Rolfs), ^He +^ He (P. Parker), e~ + ^Be (J. Bahcall), p + ^Be 
(E. Adelberger), and CNO (H. Robertson). Initial drafts of each of the sections in this paper 
were written by the discussion leaders and their close collaborators. Successive iterations 
of the paper were posted on the Internet so that they could be read and commented on by 
each member of the collaboration, resulting in an almost infinite number of iterations. Each 
section of the paper was reviewed extensively and critically by co-authors who did not draft 
that section, and, in a few cases, vetted by outside experts. 



C. Contents 

The organization of this paper reflects the organization of our workshop. Section II 
describes the theoretical justification and the phenomenological situation regarding extrap- 
olations from higher laboratory energies to lower solar energies, as well as the effects of 
electron screening on laboratory and solar fusion rates. Sections [m]-[rX| contain detailed de- 
scriptions of the current situation with regard to the most important solar fusion reactions. 
We do not consider explicitly in this review the reactions ^H(p, 7)^He, ^Li(p, a)^He, and 
®B(/5"'"z/e)*Be, which occur in the pp chain but whose rates are so fast that the precise cross 
section or decay time does not affect the energy generation or the neutrino flux calculations. 
We concentrate our discussion on those reactions that are most important for calculating 
solar neutrino fluxes or energy production. 

In our discussions at the workshop, and in the many iterations that have followed over 
the subsequent months, we placed as much emphasis on determining reliable error estimates 
as on specifying the best values. We recognize that, for applications to astronomy and to 



^ The workshop was proposed by John Bahcall, the principal editor of this paper, in a letter 
submitted to the Advisory Committee of the Institute for Nuclear Theory, August 20, 1996. W. 
Haxton, P. Parker, and H. Robertson served as joint organizers (with Bahcall) of the workshop 
and as co-editors of this paper. All of the co-authors participated actively in some stage of the 
work and/or the writing of this paper. We attempted to be complete in our review of the literature 
prior to the workshop meeting and have taken account of the most relevant work that has been 
published prior to the submission of this paper in September, 1997. 



9 



neutrino physics, it is as important to know the hmits of our knowledge as it is to record the 
preferred cross-section factors. Wherever possible, experimental results are given with la 
error bars (unless specifically noted otherwise). For a few quantities, we have also quoted 
estimates of a less precisely defined quantity that we refer to as an "effective 3a" error (or a 
maximum likely uncertainty). In order to meet the challenges and opportunities provided by 
increasingly precise solar neutrino and helioseismological data, we have emphasized in each 
of the sections on individual reactions the most important measurements and calculations 
to be made in the future. 

The sections on individual reactions, |T^-|X|, answer the questions: "What?", "How 
Well?", and "What Next?". Table | summarizes the answers to the questions "What?" and 
"How Well?" ; this table gives the best estimates and uncertainties for each of the principal 
solar fusion reactions that are discussed in greater detail later in this paper. The different 
answers to the question "What Next?" are given in the individual Sees. |l|-0. 



II. EXTRAPOLATION AND SCREENING 



A. Phenomenological Extrapolation 

Nuclear fusion reactions occur via a short-range (less than or comparable to a few fm) 
strong interaction. However, at the low energies typical of solar fusion reactions (~ 5 keV to 
30 keV), the two nuclei must overcome a sizeable barrier provided by the long-range Coulomb 
repulsion before they can come close enough to fuse. Therefore, the energy dependence of 
a (nonresonant) fusion cross section is conveniently written in terms of an S'-factor which is 
defined by the following relation: 

a(E) = ^exp{-2vrr/(E)}, (7) 

where 

m - ^ (8) 

is the Sommerfeld parameter. Here, E is the center-of-mass energy; v = (2E//i)i/2 is the 
relative velocity in the entrance channel; Zi and Z2 are the charge numbers of the colliding 
nuclei; /i = mAiA2/ {Ai + A2) is the reduced mass of the system; m is the atomic mass unit; 
and Ai and A2 are the masses (in units of m) of the reacting nuclei. 

The exponential in Eq. (^ (the Gamow penetration factor) takes into account quantum- 
mechanical tunneling through the Coulomb barrier; the exponential describes well the rapid 
decrease of the cross section with decreasing energy. The Gamow penetration factor dom- 
inates the energy dependence, derived in the WKB approximation, of the cross section in 
the low-energy limit. In the low-energy regime in which the WKB approximation is valid, 
the function S{E) is slowly varying (except for resonances) and may be approximated by 

S{E)^S{0) + S'{0)E + ^S"{0)E\ (9) 

The coefficients in Eq. (P) can often be determined by fitting a quadratic formula to labora- 
tory measurements or theoretical calculations of the cross section made at energies of order 
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100 keV to several MeV. The cross section is then extrapolated to energies, C(lOkeV), typ- 
ical of solar reactions, through Eq. (|^). However, special care has to be exercised for certain 
reactions like ''Be(p,7)^B, where the S'-factor at very low energies expected from theoretical 
considerations cannot be seen in available data (cf. discussion in Sec. |V111| ). 

The WKB approximation for the Gamow penetration factor is valid if the argument 
of the exponential is large, 2711] ^ 1. This condition is satisfied for the energies over which 
laboratory data on solar fusion reactions are usually fitted. Because the WKB approximation 
becomes increasingly accurate at lower energies, the standard extrapolation to solar-fusion 
energies is valid. 

The most compelling evidence for the validity of the approximations of Eqs. (0)- 
(P) is empirical: they successfully fit low-energy laboratory data. For example, for the 
'^He(^He,2p)'^He reaction, a quadratic polynomial fit (with only a small linear and even 
smaller quadratic term) for S{E) provides an excellent fit to the measured cross section 
over two decades in energy in which the measured cross section varies by over ten orders of 
magnitude (see discussion in Sec. 

The approximation of S{E) by the lowest terms in a Taylor expansion is supported theo- 
retically by explicit calculations for a wide variety of reasonable nuclear potentials, for which 
S{E) is found to be well approximated by a quadratic energy dependence. The specific form 
of Eq. (^ describes s-wave tunneling through the Coulomb barrier of two point-like nu- 
clei. Several well-known and thoroughly investigated effects introduce slowly-varying energy 
dependences that are not included explicitly in the standard definition of the low-energy 
S'-factor. These effects include (see, for example, Barnes, Koonin, and Langanke, 1993; 
Descouvemont, 1993; Langanke and Barnes, 1996) 1) the finite size of the colliding nuclei, 
2) nuclear structure and strong interaction effects, 3) antisymmetrization effects, 4) contri- 
butions from other partial waves, 5) screening by atomic electrons, and 6) final-state phase 
space. These effects introduce energy dependences in the S'-factor that, in the absence of 
near-threshold resonances, are much weaker than the dominant energy dependence repre- 
sented by the Gamow penetration factor. The standard picture of an S'-factor with a weak 
energy dependence has been found to be valid for the cross-section data of all nuclear re- 
actions important for the solar pp-chains. Theoretical energy dependences that take into 
account all the effects listed above are available (and have been used) for extrapolating data 
for all the important reactions in solar hydrogen burning. 

One can reduce {but not eliminate) the energy dependence of the extrapolated quantity 
by removing nuclear finite-size effects (item 1) from the data. The resulting modified S{E) 
factor is still energy dependent (because of items 2-6) and cannot be treated as a constant 
[as assumed by Dar and Shaviv (1996)]. 



B. Laboratory Screening 

It has generally been believed that the uncertainty in the extrapolated nuclear cross 
sections is reduced by steadily lowering the energies at which data can be taken in the 
laboratory. However, this strategy has some complications (Assenbaum, Langanke, and 
Rolfs, 1987) since at very low energies the experimentally measured cross section does not 
represent the bare nucleus cross section: the laboratory cross section is increased by the 
screening effects arising from the electrons present in the target (and in the projectile). The 
resulting enhancement of the measured cross section, (Texp(-E), relative to the cross section 
for bare nuclei, cr{E), can be written as 
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(10) 



Since the electron screening energy, [/g, is much smaller than the scattering energies, E, 
currently accessible in experiments, one finds (Assenbaum, Langanke, and Rolfs, 1987) 

/(E)^exp|7rr/(E)^} . (11) 

In nuclear astrophysics, one starts with the bare nuclei cross sections and corrects them 
for the screening appropriate for the astrophysical scenario (plasma screening, see Sec. |11 C| ) . 
In the laboratory experiments, the electrons are bound to the nucleus, while in the stellar 
plasma they occupy (mainly) continuum states. Therefore, the physical processes underlying 
screening effects are different in the laboratory and in the plasma. 

The enhancement of laboratory cross sections due to electron screening is well established, 
with the ^He((i, p)^He reaction being the best studied and most convincing example (Engstler 
et al, 1988; Prati et ai, 1994). However, it appeared for some time that the observed 
enhancement was larger than the one predicted by theory. This discrepancy has recently 
been removed after improved energy loss data became available for low-energy deuteron 
projectiles in helium gas. To a good approximation, atomic-target data can be corrected 
for electron screening effects within the adiabatic limit (Shoppa et ai, 1993) in which the 
screening energy, Ug, is simply given by the difference in electronic binding energy of the 
united atom and the sum of the projectile and target atoms. It appears now as if the electron 
screening effects for atomic targets can be modeled reasonably well (Langanke et ai, 1996; 
Bang et al, 1996; but see also Junker et al. 1997 ). This conclusion must be demonstrated 
for molecular and solid taigets. Experimental work on electron screening with molecular and 
solid targets is discussed in Engstler et al. (1992), while the first theoretical approaches are 
presented in Shoppa et al. (1996) (molecular) and in Boudouma, Chami, and Beaumevieille 
(1997) (sohd targets). 

Electron screening effects, estimated in the adiabatic limit, are relatively small in the 
measured cross sections for most solar reactions, including the important ^He(a,7)^Be 
and ''Be(p, 7)^B reactions (Langanke, 1995). However, both the ^He('^He, 2p)^He and the 
^^N(p, 7)^^0 data, which extend to very low energies, are enhanced due to electron screening 
and have been corrected for these effects (see Sec. ffV] and [TX]). 



C. Stellar Screening 

As shown by Salpeter (1954), the decreased electrostatic repulsion between reacting ions 
caused by the Debye-Hiickel screening leads to an increase in reaction rates. The reaction 
rate enhancement factor for solar fusion reactions is, to an excellent approximation (Gruzinov 
and Bahcall, 1998), 

/ = exp iTFT^ ' (12) 



kTR 



D 



where Ro is the Debye radius and T is the temperature. The Debye radius is defined by 
the equation Rd = {Aime'^C'^ /kT)~^/'^, where n is the baryon number density {p/rrig,^^), 

( = + (^'j^ TjiXi^i , Xi, Zi, and Ai are, respectively, the mass fraction, the 
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nuclear charge, and the atomic weight of ions of type i. The quantity /'// ^ 0.92 accounts 
for electron degeneracy (Salpeter 1954). Equation ([T^ ) is valid in the weak- screening limit 
which is defined by kTR£, ^ Z\Z2e^ ■ In the solar case, screening is weak for ZxZ2 of the 
order 10 or less (Gruzinov and Bahcall, 1998). Thus, plasma screening corrections to all 
important thermonuclear reaction rates are known with uncertainties of the order of a few 
percent. Although originally derived for thermonuclear reactions, the Salpeter formula also 
describes screening effects on the ''Be electron capture rate with an accuracy better than 1% 
(Gruzinov and Bahcall, 1997) (for ^Be(e, i^)'^Li, we have Zi = —1, and = 4). 

Two papers questioning the validity of the Salpeter formula in the weak-screening limit 
appeared during the last decade, but subsequent work demonstrated that the Salpeter for- 
mula was correct. The "3/2" controversy introduced by Shaviv and Shaviv (1996) was 
resolved by Briiggen and Gough (1997); a "dynamic screening" effect discussed by Carraro, 
Schafer, and Koonin (1988) was shown to be not present by Brown and Sawyer (1997a) and 
Gruzinov (1997). 

Corrections of the order of a few percent to the Salpeter formula come from the non- 
linearity of the Debye screening and from the electron degeneracy. There are two ways to 
treat these effects - numerical simulations (Johnson et al, 1992) and illustrative approxima- 
tions (Dzitko et al, 1995; Turck-Chieze and Lopes, 1993). Fortunately, the asymmetry of 
fluctuations is not important (Gruzinov and Bahcall, 1997), and numerical simulations of a 
spherically symmetrical approximation are possible even with nonlinear and degeneracy ef- 
fects included (Johnson et al, 1992). The discussion of intermediate screening by Graboske 
et al. (1973) is not applicable to solar fusion reactions because Graboske et al. assume 
complete electron degeneracy (cf. Dzitko et al., 1995). 

A fully analytical treatment of nonlinear and degeneracy effects is not available, but 
Brown and Sawyer (1997a) have recently reproduced the Salpeter formula by diagram sum- 
mations. It would be interesting to evaluate higher order terms (describing deviations from 
the Salpeter formula) using these or similar methods. 



III. THE pp AND pep REACTIONS 

The rates for most stellar nuclear reactions are inferred by extrapolating measurements 
at higher energies to stellar reaction energies. However, the rate for the fundamental p+p — >■ 
'^D + e'^ + i>f. reaction is too small to be measured in the laboratory. Instead, the cross section 
for the p-p reaction must be calculated from standard weak interaction theory. 

The most recent calculation was performed by Kamionkowski and Bahcall (1994), who 
used improved data on proton-proton scattering and included the effects of vacuum polar- 
ization in a self-consistent fashion. They also isolated and evaluated the uncertainties due 
to experimental errors and theoretical evaluations. 

The calculation of the p-p rate requires the evaluation of three main quantities: (i) the 
weak- interaction matrix element, (ii) the overlap of the pp and deuteron wave functions, and 
(iii) mesonic exchange-current corrections to the lowest-order axial- vector matrix element. 

The best estimate for the logarithmic derivative, 

S'{0) = (11.2 ± 0.1) MeV-\ (13) 

is still that of Bahcall and May (1968). At the Gamow peak for the pp reaction in the 
Sun, this linear term provides only an 0{1%) correction to the E = value. The quadratic 
correction is several orders of magnitude smaller, and therefore negligible. Furthermore, 
the 1% uncertainty in Eq. ( [T3| ) gives rise to an C(0.01%) uncertainty in the total reaction 
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rate. This is negligible compared with the uncertainties described below. Therefore, in the 
following, we focus on the E = cross-section factor. 

At zero relative energy, the S'-factor for the pp reaction rate can be written (Bahcall and 
May, 1968, 1969) , 

^(0) = 6vr^m,caln2 ^ (^)' (1 + 5)^ (14) 

where a is the fine-structure constant; rUp is the proton mass; Gy and Ga are the usual Fermi 
and axial-vector weak coupling constants; 7 = {2^EdY^'^ = 0.23161 fm~^ is the deuteron 
binding wave number (/i is the proton-neutron reduced mass and is the deuteron binding 
energy); is the phase-space factor for the pp reaction (Bahcall, 1966) with radiative 
corrections; (/t)o+^o+ is the ft value for superallowed 0"^ 0+ transitions (Savard et al, 
1995); A is proportional to the overlap of the pp and deuteron wave functions in the impulse 
approximation (to be discussed below); and 6 takes into account mesonic corrections. 
Inserting the current best values, we find 

S(0) = 4.00 X 10- MeV b (W^) " (J^) (^) ' (JL) (111] ^ (ig) 

V 3073 sec ; V6.92; 1.2654 / 1^0.144/ 1.01 / ^ ^ 

We now discuss the best estimates and the uncertainties for each of the factors which appear 
in Eq. (|15D. 

The quantity A^ is proportional to the overlap of the initial-state pp wave function and 
the final-state deuteron wave function. The wave functions are determined by integrating the 
Schrodinger equations for the two-nucleon systems with an assumed nuclear potential. The 
two-nucleon potentials cannot be determined from first principles, but the parameters in any 
given functional form for the potentials must fit the experimental data on the two-nucleon 
system. By trying a variety of dramatically different functional forms, we can evaluate the 
theoretical uncertainty in the final result due to ignorance of the form of the two-nucleon 
interaction. 

The proton-proton wave function is obtained by solving the Schrodinger equation for 
two protons that interact via a Coulomb plus nuclear potential. The potential must fit 
the pp scattering length and effective range determined from low-energy pp scattering. In 
Kamionkowski and Bahcall (1994), five forms for the nuclear potential were considered: a 
square well, Gaussian, exponential, Yukawa, and a repulsive-core potential. The uncertainty 
in A^ from the pp wave function is small because there is only a small contribution to the 
overlap integral from radii less than a few fm (where the shape of the nuclear potential 
affects the wave function). At larger radii, the wave function is determined by the measured 
scattering length and effective range. The experimental errors in the pp scattering length 
and effective range are negligible compared with the theoretical uncertainties. 

Similarly, the deuteron wave function must yield calculated quantities consistent with 
measurements of the static deuteron parameters, especially the binding energy, effective 
range, and the asymptotic ratio of D- to S'-state deuteron wave functions. In Kamionkowski 
and Bahcall (1994), seven deuteron wave functions which appear in the literature were 
considered. The spread in A due to the spread in assumed neutron-proton interactions was 
0.5%, and the uncertainty due to experimental error in the input parameters was negligible. 

Figure |l| shows why the details of the nuclear physics are unimportant. The figure 
displays the product of the radial pp and deuteron wave functions, Upp{r) and Ud{r). The 
wavelength of the pp system is more than an order of magnitude larger than the extent of 
the deuteron wave function, so the shape of the curve shown in Fig. |l] is independent of pp 
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energy. Most of the contribution to the overlap integral between the pp wave function and the 
deuteron wave function comes from relatively large radii where experimental measurements 
constrain the wave function most strongly. The assumed shape of the nuclear potential 
produces visible differences in the wave function only for r ^ 5 fm, and these differences are 
small. Furthermore, only ~ 40% of the integrand comes from r ^ 5 fm and ~ 2% of the 
integrand comes from r ^ 2 fm. 

Including the effects of vacuum polarization and the best available experimental param- 
eters for the deuteron and low-energy pp scattering, one finds (Kamionkowski and Bahcall, 



where the first uncertainty is due to experimental errors, and the second is due to theoretical 
uncertainties in the form of the nuclear potential. 

An anomalously high value of = 7.39 was obtained by Gould and Guessoum (1990), 
who did not make clear what values for the pp scattering length and effective range they 
used. Even by surveying a wide variety of nuclear potentials that fit the observed low-energy 
pp data, Kamionkowski and Bahcall (1994) never found a value of A^ greater than 7.00. We 
therefore conclude that the large value of A^ reported by Gould- Guessoum is caused by 
either a numerical error or by using input data that contradict the existing pp scattering 
data. 

The calculation of A^ includes the overlap only of the s-wave (i.e., orbital angular mo- 
mentum / = 0) part of the pp wave function and the S state of the deuteron. Because the 
matrix element is evaluated in the usual allowed approximation, D— state components in 
the deuteron wave function do not contribute to the transition. 

We use (/t)o+^o+ = (3073.1 ± 3.1), which is the ft value for superallowed 0+ 0"*" 
transitions that is determined from experimental rates corrected for radiative and Coulomb 
effects (Savard et al, 1995). This value is obtained from a comprehensive analysis of data 
on numerous O"*" 0^ superallowed decays. After radiative corrections, the ft values for all 
such decays are found to be consistent within the quoted error. 

Barnett et al. (1996) recommend a value Ga/Gv = 1.2601 ±0.0025, which is a weighted 
average over several experiments that determine this quantity from the neutron decay asym- 
metry. However, a recent experiment (Abele et al., 1997) has obtained a slightly higher value. 
We estimate that if we add this new result to the compilation of Barnett et al. (1996), the 
weighted average will be Ga/Gv = 1.2626 ±0.0033. Alternatively, Ga/Gv may be obtained 
from (/t)o+^o+ and the neutron /t- value from 



For the neutron lifetime, we use tn = (888 ± 3) sec. The range spanned by this central value 
and the la uncertainty covers the ranges given by the recommended value and uncertainty 
(887 ± 2.0) of Barnett et al. (1996) and the value and uncertainty (889.2 ± 2.2), obtained if 
the results of Mampe (1993) — which have been called into question by Ignatovich (1995) — 
are left out of the compilation. We use the neutron phase-space factor, /„ = 1.71465 
(including radiative corrections), obtained in Wilkinson (1982). Inserting the ft values into 
Eq. (|T3), we find Ga/Gv = 1.2681 ± 0.0033, which is slightly larger (by 0.0055 or 0.4%) 
than the value obtained from neutron decay distributions. To be conservative, we take 
Ga/Gv = 1.2654 ± 0.0042. 

Considerable work has been done on corrections to the nuclear matrix element for the 
exchange of vr and p mesons (Gari and Huffman, 1972; Dautry, Rho, and Riska, 1976), which 



1994) 



A^ = 6.92 X (1 ±0.002+°;°^^), 



(16) 




(17) 
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arise from nonconscrvation of the axial- vector current. By fitting an effective interaction 
Lagrangian to data from tritium decay, one can show phenomenologically that the mesonic 
corrections to the pp reaction rate should be small (of order a few percent) (Blin-Stoyle 
and Papageorgiou, 1965). Heuristically, this is because most of the overlap integral comes 
from proton-proton separations that are large compared with the typical (~ 1 fm) range 
of the strong interactions. In tritium decay, most of the overlap of the initial and final 
wave functions comes from a much smaller radius. If mesonic effects are to be taken into 
account properly, they must be included self-consistently in the nuclear potentials inferred 
from data and in the calculation of the overlap integral described above. Here, we advocate 
following the conservative recommendation of Bahcall and Pinsonneault (1992) in adopting 
5 — O.Olloioi- The central value is consistent with the best estimates from two recent 
calculations which take into account p as well as tt exchange (Bargholtz, 1979; Carlson et 
al, 1991) . 

The quoted error range for 5 could probably be reduced by further investigations. The 
primary uncertainty is not in the evaluation of exchange current matrix elements, since the 
deuteron wave function is well determined from microscopic calculations, but in the meson- 
nucleon-delta couplings that govern the strongest exchange currents. The coupling constant 
combinations appearing in the present case are similar to those contributing to tritium beta- 
decay, another system for which accurate microscopic calculations can be made. Thus the 
mcasiu^cd '^H lifetime places an important constraint on the exchange current contribution 
to the pp reaction. In the absence of a detailed analysis of this point, the error adopted 
above, which spans the range of recently published calculations, remains appropriate. But 
we point out that the '^H lifetime should be exploited to reduce this uncertainty. 

For the phase-space factor f^, we have taken the value without radiative corrections, 
fpp = 0.142 (Bahcall and May, 1969) and increased it by 1.4% to take into account ra- 
diative corrections to the cross section. Although first-principle radiative corrections for 
this reaction have not been performed, our best ansatz (Bahcall and May, 1968) is that 
they should be comparable in magnitude to those for neutron decay (Wilkinson, 1982). To 
obtain the magnitude of the correction for neutron decay, we simply compare the result 
= 1.71465 with radiative corrections obtained in Wilkinson (1982) to that obtained 
without radiative corrections in Bahcall (1966). We estimate that the total theoretical un- 
certainty in this approximation for the pp phase-space factor is 0.5%. Therefore, we adopt 
fpp — 0.144 X (1 ± 0.005), where the error is a total theoretical uncertainty (see Bahcall, 
1989). It would be useful to have a first-principles calculation of the radiative corrections 
for the pp interaction. 

Amalgamating all these results, we find that the current best estimate for the pp cross- 
section factor, taking account of the most recent experimental and theoretical data, is 

5(0) = 4.00 X 10"^^ (1 ± Q.mtQ^) MeV b, (18) 

where the first uncertainty is a la experimental error, and the second uncertainty is one-third 
the estimated total theoretical uncertainty. 

Ivanov et al. (1997) have recently calculated the pp reaction rate using a relativistic field 
theoretic model for the deuteron. Their calculation is invalidated by, among other things, 
the fact that they used a zero-range effective interaction for the protons, in confiict with 
low-energy pp scattering experiments (see Bahcall and Kamionkowski, 1997). 

The rate for the p -|- e" — > -|- i/g reaction is proportional to that for the pp reaction. 
Bahcall and May (1969) found that the pep rate could be written, 

Rp^p ~ 5.51 X 10-^(1 + X)T^^'^{1 + QmTQ)Rpp, (19) 
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where p is the density in g cm~^, X is the mass fraction of hydrogen, Tg is the temperature 
in units of 10^ K, and Rpp is the pp reaction rate. This approximation is accurate to 
approximately 1% for the temperature range, 10 < < 16, relevant for solar-neutrino 
production. Therefore, the largest uncertainty in the pep rate comes from the uncertainty 
in the pp rate. 



IV. THE ='He(='He, 2p)'^U.e REACTION 

The solar Gamow energy of the 3He(3He, 2p)^He reaction is at Eq = 21.4 keV (see Eq. 1). 
As early as 1972, there were desperate proposals (Fetisov and Kopysov, 1972; Fowler, 1972) 
to solve the solar neutrino problem^ that suggested a narrow resonance may exist in this 
reaction at low energies. Such a resonance would enhance the ^He + ^He rate at the expense 
of the ^He + '^He chain, with important repercussions for production of ^Be and neutrinos. 
Many experimental investigations [see Rolfs and Rodney (1988) for a list of references] have 
searched for, but not found, an excited state in ^Be at £^2; ~ 11.6 MeV that would correspond 
to a low-energy resonance in ^He + ^He. Microscopic theoretical models (Descouvemont, 
1994; Csoto, 1994) have also shown no sign of such a resonance. 

Microscopic calculations of the ^He('^He,2p)^He reaction (Vasilevskii and Rybkin, 1989; 
Typel et al, 1991) view this reaction as a two-step process: After formation of the com- 
pound nucleus, the system decays into an a-particle and a 2-proton cluster. The latter, 
being energetically unbound, finally decays into two protons. This, however, is expected 
to occur outside the range of the nuclear forces. In Typel et al. (1991), the model space 
was spanned by ^He-|-2p and ^He-|-^He cluster functions as well as configurations involving 
■^He pseudostates. The calculation reproduces the measured S-factors for E < 300 keV 
reasonably well and predicts S{0) ~ 5.3 MeV b, in agreement with the measurements dis- 
cussed later in this section. Further confidence in the calculated energy dependence of the 
low-energy ^He(^He,2p)^He cross sections is gained from a simultaneous microscopic cal- 
culation of the analog '^H(^H, 2n)^He reaction, which again reproduces well the measured 
energy dependence of the ^H-|-^H fusion cross sections (Typel et al, 1991). Recently, De- 
scouvemont (1994) and Csoto (1997b, 1998) have extended the microscopic calculations to 
include ^Li -|- p configurations. Their calculated energy dependences, however, are in slight 
disagreement with the data, possibly indicating the need for a genuine 3-body treatment of 
the final continuum states. 

The relevant cross sections for the ^He(^He, 2p)^B.e reaction have recently been measured 
at the energies covering the Gamow peak. The data have to be corrected for laboratory 
electron screening effects. Note that the extrapolation given by Krauss et al. (1987) and 
used in Dar and Shaviv (1996) (5'(0) = 5.6 keV b) is too high, because it is based on 
low-energy data that were not corrected for electron screening. 

The reaction data show that at energies below 1 MeV the reaction proceeds predomi- 
nately via a direct mechanism and that the angular distributions approach isotropy with 
decreasing energy. The energy dependence of <j{E) — or equivalently of the cross-section 
factor S{E) — observed by various groups (Bacher and Tombrello, 1965; Wang et al, 1966; 
Dwarakanath and Winkler, 1971; Dwarakanath, 1974; Krauss, Becker, Trautvetter, and 



^In 1972, the "solar neutrino problem" consisted entirely of the discrepancy between the predicted 
and measured rates in the Homestake experiment (see Bahcall and Davis, 1976). 



17 



Rolfs, 1987; Greife et al, 1994; Arpesella et al, 1996; Junker et al, 1997) presents a con- 
sistent picture. The only exception is the experiment of Good, Kunz, and Moak (1951), for 
which the discrepancy is most likely caused by target problems (^He trapped in an Al foil). 

The absolute S{E) values of Dwarakanath and Winkler (1971), Krauss, Becker, Trautvet- 
ter, and Rolfs (1987), Greife et al. (1994), Arpesella et al. (1996), and Junker et al. (1997) 
are in reasonable agreement, although they are perhaps more consistent with a systematic 
uncertainty of 0.5 MeV b. The data of Wang et al. (1966) and Dwarakanath (1974) are 
lower by about 25%, suggesting a renormalization of their absolute scales. However, in view 
of the relatively few data points reported in Wang et al. (1966) and Dwarakanath (1974), 
and their relatively large uncertainties — in comparison to other data sets — we suggest that 
the data of Wang el al. (1966) and Dwarakanath (1974) can be omitted without significant 
loss of information. 

Figure ^ is adapted from Fig. 9 of Junker et al. (1997). The data shown are from 
Dwarakanath and Winkler (1971), Krauss, Becker, Trautvetter, and Rolfs (1987), Arpesella 
et al. (1996), and Junker et al. (1997) . The data provide no evidence for a hypothetical 
low-energy resonance over the entire energy range that has been investigated experimentally. 

Because of the effects of laboratory atomic electron screening (Assenbaum, Langanke, 
and Rolfs, 1987), the low-energy '^He(^He, 2p)^He measurements must be corrected in order 
to determine the "bare" nuclear S'-factor. Assume, for specificity, a constant laboratory 
screening energy of Ue = 240 eV, corresponding to the adiabatic limit for a neutral ^He beam 
incident on the atomic ^He target. If we assume that the projectiles are singly ionized, the 
adiabatic screening energy increases only slightly to Ue ~ 250 eV. TDHF calculations for 
atomic screening of low-Z targets (Shoppa et al, 1993; Shoppa et al, 1996) have shown that 
the adiabatic limit is expected to hold well at the low energies where screening is important. 
Junker et al. (1997) have converted published laboratory measurements 5'iab(-E) to bare 
nuclear S-factors S{E) using the relation S{E) = S'iab(-E) exp{—7!-ri{E)Ue/E), with Ue = 240 
eV [cf. Eq. (|10|) and Eq. {^]. 

The resulting bare ^-factors were fit to Eq. (^). Junker et al. (1997) find ^(O) = 5.40 ± 
0.05 MeV b, S'{0) = -4.1 ± 0.5 b, and S'' (0) = 4.6 ± 1.0 b/MeV, but important systematic 
uncertainties must also be included as in Eq. (|20D below. An effective 3a uncertainty of 
about ±0.30 MeV b due to lack of understanding of electron screening in the laboratory 
experiments should be included in the error budget for 5(0) (cf. Junker et al, 1997). 

The cross-section factor at solar energies is relatively well known by direct measurements 
(see Fig.0). Junker et al. (1997) give 

S{Eo) = 5.3 ± 0.05(stat) ± 0.30(syst) ± 0.30(Ue) MeV b, (20) 

where the first two quoted la errors are from statistical and systematic experimental un- 
certainties and the last error represents a maximum likely error (or effective 3a error) due 
to the lack of complete understanding of laboratory electron screening. The data seem to 
suggest that the effective value of Ue may be larger than the adiabatic limit. 

Future experimental efforts should extend the S{E) data to energies at the low-energy 
tail of the solar Gamow peak, i.e. at least as low as 15 keV. Furthermore, improved data 
should be obtained at energies from E = 25 keV to 60 keV to confirm or reject the possibility 
of a relatively large systematic error in the S{E) data near these energies. On the theoretical 
side, an improved microscopic treatment is highly desirable. 
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V. THE =^He(Q;,7)''Be REACTION 



The relative rates of the ^He(a,7)^Be and '^He(^He, 2p)'^He reactions determine what 
fractions of pp-chain terminations result in ^Be or neutrinos. 

Since the ^He(Q!,7)^Be reaction at low energies is essentially an external direct-capture 
process (Christy and Duck, 1961), it is not surprising that direct-capture model calculations 
of different sophistication yield nearly identical energy dependences of the low-energy S- 
factor. Both the microscopic cluster model (Kajino and Arima, 1984) and the microscopic 
potential model (Langanke, 1986) correctly predicted the energy dependence of the low- 
energy '^H(a, 7)'^Li cross section (the isospin mirror of ^B.e{a, 7)^Be ) before it was precisely 
measured by Brune, Kavanagh, and Rolfs (1994). The absolute value of the cross section 
was also predicted to an accuracy of better than 10% from potential model calculations by 
Langanke (1986) and Mohr et al. (1993). 

Separate evaluations of this energy dependence based on the Resonating Group Method 
(Kajino, Toki, and Austin, 1987) and on a direct-capture cluster model (Tombrello and 
Parker, 1963) agree to within ±1.25% and are also in good agreement with the measured 
energy dependence (see also Igamov, Tursunmuratov, and Yarmukhamedov, 1997). This 
confluence of experiments and theory is illustrated in Fig. ^. Even more detailed calculations 
are now possible (cf. Csoto, 1997a). 

Thus the energy dependence of the '^He(Q;, 7)^Be reaction seems to be well determined. 
The only free parameter in the extrapolation to thermal energies is the normalization of 
the energy dependence of the cross sections to the measured data sets. While the energy 
dependence predicted by the existing theoretical models is in good agreement with the 
energy dependence of the measured cross sections, it would be useful to explore how robust 
this energy dependence is for a wider range of models. Extrapolations based on physical 
models should be used; such extrapolations are more credible than those based only on the 
extension of multiparameter mathematical fits (e.g., those of Castellani et al., 1997). 

There are six sets of measurements of the cross section for the ^He(Q;,7)^Be reaction 
that are based on detecting the capture gamma rays (Table 0). The weighted average of 
the six prompt 7-ray experiments yields a value of 5*34(0) = (0.507 ± .016) keV b, based on 
extrapolations made using the calculated energy dependence for this direct-capture reaction. 
In computing this weighted average, we have used the renormalization of the Krawinkel et 
al (1982) result by Hilgemeier et al (1988). 

There are also three sets of cross sections for this reaction that are based on measurements 
of the activity of the synthesized ''Be (Table y). These decay measurements have the 
advantage of determining the total cross section directly, but have the disadvantage that 
(since the source of the residual activity can not be uniquely identified) there is always the 
possibility that some of the ^Be may have been produced in a contaminant reaction that 
evaded background tests. The three activity measurements (when extrapolated in the same 
way as the direct-capture gamma ray measurements) yield a value of S'34(0) = (0.572±0.026) 
keV b, which differs by about 2.5(T with the value based on the direct-capture gamma rays. 

It has been suggested that the systematic discrepancy between these two data sets might 
arise from a small monopole (EO) contribution to which the prompt measurements would be 
much less sensitive and whose contribution could have been overlooked. However, estimates 
of the EO contribution are consistently found to be exceedingly small in realistic models of 
this reaction; they are of order a^, whereas the leading contribution is of order a (the fine 
structure constant). The importance of any EO contributions would be further suppressed 
by the fact that they would have to come from the p-wave incident channel, in contrast to 
the s-wave incident channel which is responsible for the dominant El contributions. (See 
Fig. |.) 
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When the nine experiments are combined, the weighted mean is 6*34 (0) = (0.533 ± 
0.013) keV b, with = 13.4 for 8 degrees of freedom. The probabihty of such a distribution 
arising by chance is 10%, and that, together with the apparent grouping of the results 
according to whether they have been obtained from activation or prompt-gamma yields, 
suggests the possible presence of a systematic error in one or both of the techniques. An 
approach that gives a somewhat more conservative evaluation of the uncertainty is to form 
the weighted means within each of the two groups of data (the data show no indication of 
non-statistical behavior within the groups), and then determine the weighted mean of those 
two results. In the absence of information about the source and magnitude of the excess 
systematic error, if any, an arbitrary but standard prescription can be adopted in which the 
uncertainties of the means of the two groups (and hence the overall mean) are increased by 
a common factor of 3.7 (in this case) to make = 0.46 for 1 degree of freedom, equivalent 
to making the estimator of the weighted population variance equal to the weighted sample 
variance. The uncertainty in the extrapolation is common to all the experiments, although 
it is likely to be only a relatively minor contribution to the overall uncertainty. The result 
is our recommended value for an overall weighted mean: 

^34(0) = 0.53 ± 0.05 keV b. (21) 

The slope, >S"(0), is well determined within the accuracy of the theoretical calculations (e.g., 
Parker and Rolfs, 1991): 

S'{0) = -0.00030 b. (22) 

Neither the theoretical calculations nor the experimental data are sufficiently accurate to 
determine a second derivative. 

Dar and Shaviv (1996) quote a value of -S'34(0) — 0.45 keV b, about 1.5(7 lower than our 
best estimate. The difference between their value and our value can be traced to the fact 
that Dar and Shaviv fit the entire world set of data points as a single group to obtain one 
5*34(0) intercept, rather than fitting independently each of the nine independent experiments 
and then combining their intercepts to determine a weighted average for <S'34(0). The Dar 
and Shaviv method thereby overweights the experiments of Krawinkel et al. (1982) and 
Parker and Kavanagh (1963) because they have by far the largest number of data points 
(39 and 40, respectively) and underweights those experiments which have only 1 or 2 data 
points (e.g., the activity measurements). Systematic uncertainties, such as normalization 
errors, common to all the points in one data set make it impossible to treat the common 
points as statistically independent and uncorrelated, and thus the Dar and Shaviv method 
distorts the average. 



VI. THE ^He(p, e+ -|- i/e)^He REACTION 

The hep reaction, 

^He + p ^ ^He + e+ + u^, (23) 

produces neutrinos with an endpoint energy of 18.8 MeV, the highest energy expected for 
solar neutrinos. The region between 15 MeV and 19 MeV, above the endpoint energy for 
neutrinos and below the endpoint energy for hep neutrinos, is potentially useful for solar 
neutrino studies since the background in electronic detectors is expected to be small in this 
energy range. For a given solar model, the flux of hep neutrinos can be calculated accurately 
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once the S*— factor for reaction (^) is specified. The rate of the hep reaction is so slow that 
it does not affect the solar structure calculations. The calculated hep flux is very small 
(~ 10^ cm~^s~^, Bahcall and Pinsonneault, 1992), but the interaction cross section is so 
large that the hep neutrinos are potentially detectable in sensitive detectors hke SNO and 
Superkamiokande (Bahcall, 1989). 

The thermal neutron cross section on ^He has been measured accurately in two separate 
experiments (Wolfs et al, 1989; Wervelman, et al, 1991). The results are in good agreement 
with each other. 

Unfortunately, there is a complicated relation between the measured thermal-neutron 
cross section and the low-energy cross-section factor for the production of hep neutrinos. 
The most recent detailed calculation (Schiavilla, et al., 1992) that includes A-isobar degrees 
of freedom yields low-energy cross-section factors calculated, with specific assumptions, in 
the range ^(0) = 1.4 x lO'^o keV b to ^(0) = 3.2 x IQ-^o keV b . Less sophisticated 
calculations yield very different answers (see Wolfs et al., 1989; Wervelman et al, 1991; see 
also the detailed calculation by Carlson et al, 1991). 

There are significant cancellations among the various matrix elements of the one- and 
two-body parts of the axial current operator. The inferred S'-factor is particularly sensitive 
to the model for the axial exchange-current operator. The uncertainties in the various 
components of the exchange- current operator and the uncertainty in the weak coupling 
constant Qijna introduce a substantial uncertainty in 5'(0). Schiavilla et al. use different 
input parameters that reflect these uncertainties, and provide a range of calculated >S'(0). 

We adopt as a best estimate low-energy cross-section factor a value in the middle of the 
range calculated by Schiavilla et al. (1992), 

S{0) = 2.3 X lO'^^keV b. (24) 

There is no satisfactory way of determining a rigorous error to be associated with this 
best estimate. However, we note that a factor of 2.5, up or down, spans the entire range 
of theoretical estimates that are in the literature and therefore corresponds to the "total 
theoretical error" often used in solar neutrino studies (Bahcall, 1989) as a substitute for a 
rigorously determined 3a uncertainty. 



Theoretical studies that could predict the cross-section factor for reaction (^) with 
greater accuracy would be important since the hep neutrino flux contains significant infor- 
mation about both solar fusion and neutrino properties. 



VII. ^Be ELECTRON CAPTURE 

The ''Be electron capture rate under solar conditions has been calculated using an explicit 
picture of continuum-state and bound-state electrons and independently using a density 
matrix formulation that does not make assumptions about the nature of the quantum states. 
The two calculations are in excellent agreement within a calculational accuracy of about 1%. 

The fluxes of both ^Be and ^B solar neutrinos are proportional to the ambient density 
of ^Be ions. The flux of ^Be neutrinos, 0(''Be), depends upon the rate of electron capture, 
R{e), and the rate of proton capture, R{p), only through the ratio 

0(7Be) oc , ^^""^ , (25) 
' R{e)^R{p) ^ ' 

With standard parameters, solar models yield R{j)) ~ 10~^i?(e). Therefore, Equation (^5|) 
shows that the flux of ''Be neutrinos is actually independent of the local rates of both the 



21 



electron capture and the proton capture reactions to an accuracy of better than 1% . The 
^Be flux depends most strongly on the branching between the ^He-^He and the '^He-^He 
reactions. The ®B neutrino flux is proportional to R{p)/[R{e) + R{p)] and therefore the 
flux is inversely proportional to the electron-capture rate. 

The first detailed calculation of the ^Be electron capture rate from continuum states 
under stellar conditions was by Bahcall (1962), who considered the thermal distribution of 
the electrons, the electron-nucleus Coulomb effect, relativistic and nuclear size corrections, 
and a numerical self-consistent Hartree wave function needed for evaluating the electron 
density at the nucleus in laboratory decay (for comparison with the electron density in 
stars). Iben, Kalata, and Schwartz (1967) made the first explicit calculation of the effect of 
bound electron capture. They included the effects of the stellar plasma in the Debye-Hiickel 
approximation and demonstrated that electron screening decreases significantly the bound 
rate compared to the case where screening is neglected. 

Applying the same Debye-Hiickel screening picture to continuum states, Bahcall and 
Moeller (1969) showed that plasma effects on the continuum capture rate were small. Bah- 
call and Moeller (1969) also formulated the total capture rate in a convenient analytic form, 
which is in general use today (Bahcall, 1989), and averaged the capture rates over three dif- 
ferent solar models. Let R = R{e) be the total capture rate and C be the rate of capture from 
the continuum only. Bahcall and Moeller (1969) found that the ratio of total rate to contin- 
uum rate averaged over the solar models was < i?/C > ~ <C/R>~^ = 1.205 ± 0.005. 

Watson and Salpeter (1973) first drew attention to the small number (~ 3) of ions per 
Debye sphere in the solar interior; they emphasized the possible importance of thermal 
plasma fluctuations on the bound-state electron-capture rate. Johnson et al. (1992) per- 
formed a series of detailed calculations, especially for the bound state capture rate, using a 
form of self-consistent Hartree theory. They derived a correction to the usual total rate of 
about 1.4%. 

Using the previously-calculated electron capture rate as a function of temperature, den- 
sity, and composition, Bahcall (1994) calculated the fraction of decays from bound states 
and found that the ratio of total to continuum captures was R/C = 1.217 ± 0.002 for three 
modern solar models, which is about 1% larger than the results of Bahcall and Moeller 
(1969) cited earlier. Using this slightly higher bound-state fraction, we find 

RCBe + e-) = 5.60 x 10"9(p//ie)T6"^/^[l + 0.004(T6 - 16)]s-\ (26) 

where fie is the electron mean molecular weight. In most recent discussions (Bahcall and 
Moeller, 1969; Bahcall, 1989), the numerical coefficient in Eq. (p6|) was 5.54 instead of 5.60. 
The slightly higher value given here reflects the newer determination of the bound fraction 
(Bahcall, 1994). 

Most recently, Gruzinov and Bahcall (1997) abandoned the standard picture of bound 
and continuum states in the solar plasma and have instead calculated the total electron 
capture rate directly from the equation for the density matrix (Feynman, 1990) of the 
plasma. Their numerical results agree to within 1% with the standard result obtained with 
an explicit picture of bound and continuum electron states. They also show that a simple 
heuristic argument, derivable from the density matrix formulation, gives an analytic form 
for the effect of the solar plasma that is of the familiar Salpeter (1954) form and agrees 
to within 1% with the numerical calculations.^ An explicit Monte Carlo calculation of the 
effects of fluctuations, not required to be spherically symmetric, shows that the net effect of 



^Even more recently, Brown and Sawyer (1997b) have re-investigated the electron capture prob- 
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fluctuations is less than 1% of the total capture rate. This result is surprising given the small 
number of ions in the Debye sphere (Watson and Salpeter, 1973). However, the fact that 
fluctuations are unimportant can be understood (or at least made plausible) using second- 
order perturbation theory in the density-matrix formulation. The effect of fluctuations is 
indeed shown (Gruzinov and Bahcall, 1997) to depend upon an inverse power (—5/3) of the 
number of ions in the Debye sphere. But, the dimensionless coefficient is tiny (2 x 10~^). 
The net result of the calculations performed with the density-matrix formalism is to confirm 
to high accuracy the standard ''Be electron-capture rate given here in Eq. (pU]). 

How accurate is the present theoretical capture rate, Rl The excellent agreement be- 
tween the numerical results obtained using different physical pictures (models for bound and 
continuum states and the density matrix formulation) suggests that the theoretical capture 
rate is relatively accurate. Moreover, a simple physical argument shows (Gruzinov and Bah- 
call, 1997) that the effects of electron screening on the total capture rate can be expressed 
by a Salpeter factor (Salpeter, 1954) that yields the same numerical results as the detailed 
calculations. The simplicity of this physical argument provides supporting evidence that the 
calculated electron capture rate is robust. 

The largest recognized uncertainty arises from possible inadequacies of the Debye screen- 
ing theory. Johnson et al. (1992) have performed a careful self-consistent quantum mechan- 
ical calculation of the effects on the ^Be electron capture rate of departures from the Debye 
screening. They conclude that Debye screening describes the electron capture rates to within 
2%. Combining the results of Gruzinov and Bahcall (1997) and of Johnson et al. (1992), 
we conclude that the total fractional uncertainty, 6R/R, is small and that (at about the la 
level) 

6RCBe + e-) 
fi('Be + e-) ^ 



VIII. THE ^Be(p,7) *B REACTION 



A. Introduction 



The neutrino event rate in all existing solar neutrino detectors, except for those based 
on the ^^Ga(i/, e) reaction, is either dominated by (in the case of the Homestake Mine ^^Cl 
detector), or almost entirely due to (in the cases of the Kamiokande, Super-Kamiokande, and 
SNO detectors), the high-energy neutrinos produced in ^B decay. It is therefore important 



lem using multi-particle field-theory methods. Their technique automatically gives the correct 
weighting with Fermi statistics (a small correction) including an account of bound states which 
obviates the need for "Saha-like" reasoning. They derive analytic sum rules which confirm the 
Gruzinov and Bahcall (1997) result that the Salpeter (1954) correction holds to good accuracy in 
the electron capture process. 
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to assess critically the information needed to predict the solar production of ^B.Q The most 
poorly known quantity in the entire nucleosynthetic chain that leads to is the rate of 
the final step, the ''Be(p, 7)^8 reaction which has a Q-value of 137.5 ± 1.2 keV (Audi and 
Wapstra, 1993). 

The '^Be(p, 7)^B rate is conventionally given in terms of the zero-energy 5'-factor S'i7(0). 
This quantity is deduced by extrapolating the measured absolute cross sections, which have 
been studied to energies as low as Ep = 134 keV, to the astrophysically relevant regime. 

Due to the small binding energy of ^B, the ^Be(p, 7)^B reaction at low energies is an 
external, direct-capture process (Christy and Duck, 1961). Consequently the energy de- 
pendence of the S'-factor for E < 300 keV is almost model-independent (Williams and 
Koonin, 1981; Csoto, 1997a; Timofeyuk, Baye, and Descouvemont, 1997) and is given by 
the predicted ratio of El capture from ''Be + p s-waves and waves into the ^B ground state 
(Robertson, 1973; Barker, 1980). The S'-factor is expected to exhibit a modest rise at solar 
energies due to the energy dependences of the Whittaker asymptotics of the ground state, 
the regular Coulomb functions describing the ^Be + p scattering states, and the E^ dipole 
phase-space factor. Because this expected rise of the S'-factor towards solar energies cannot 
be seen at the energies at which capture data is currently available, extrapolations that do 
not incorporate the correct physics of the low-energy ''Be(p, 7)'^B reaction, for example, the 
extrapolation presented by Dar and Shaviv (1996), are not correct. 

We have fitted Johnson et al.'s (1992) microscopic calculations of S{E) to quadratic 
functions between 20 keV and 300 keV. The overall normalization was allowed to float 
and only the energy dependence was fitted. The results were practically the same for the 
Minnesota force (Chwieroth et ai, 1973) and the Hasegawa-Nagata force (Furutani et al, 
1980). A combined fit, weighting the results from both force laws equally, yields S'(0) / S(0) = 
-0.7 ±0.2 MeV"^ and S"{0)/2S{0) = 1.9 ±0.3 MeV"^ which are our recommended values. 
The quadratic formulae given above reproduce the detailed microscopic calculations to an 
accuracy of ±0.3 eV b in the energy range to 300 keV. 

At moderate energies, say E > 400 keV, the ''Be(p, 7)'^B S'-factor becomes model- 
dependent (e.g., Csoto, 1997a), because at these energies the capture process probes the 
internal ®B wave function and becomes sensitive to nuclear structure. The argument of 
Nunes, Crespo, and Thompson (1997) that the energy dependence of Sn is sensitive to core 
polarization effects has been found to be invalid and the paper has been withdrawn by the 
authors. At the present time, statistical and systematic errors in the experimental data 
dominate the uncertainty in the low-energy cross-section factor (see also Turck-Chieze et 
ai, 1993). A measurement of the cross section below 300 keV with an uncertainty signifi- 
cantly better than 5% would make a major contribution to our knowledge of this reaction. 
A measurement of the ^Be quadrupole moment would also help distinguish between different 
nuclear models for the ''Be(p, 7)^B reaction (see Csoto et ai, 1995). 

We begin by reviewing the history of direct measurements of the ^Be(p, 7)^B cross section. 
Then we discuss recent indirect attempts to determine the cross section. Finally we make 
recommendations for S'i7(0). 



^The shape of the energy spectrum from decay is the same (Bahcall, 1991), to one part in 10^, 
as the shape determined by laboratory experiments and is relatively well known (see Bahcall et 
al, 1996). 
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B. Direct ^Be(p, 7)^B measurements 

The first experimental study of ^Be(p, 7)^6 was made by Kavanagli (1960) wlio detected 
the (3^ activity. This pioneering measurement was followed by an experiment by Parker 
(1966, 1968), who improved the signal-to-background by detecting the /3-delayed a's, a strat- 
egy followed in all subsequent works. Subsequently, extensive measurements were reported 
by Kavanagh et al. (1969) in the energy region Ep = 0.165 to 10 MeV, and by Vaughn et al. 
(1970) at 20 proton energies between 0.953 and 3.281 MeV. The most recent published works 
are a single point at Ep = 360 keV by Wiczorck et al. (1977) and a very comprehensive 
and careful experiment by Filippone et al. (1983a,b), who measured the cross section at 25 
points at center-of-mass energies between 0.117 and 1.23 MeV. The cross section displays a 
strong J'^ = 1+ resonance at Ep = 0.72 MeV, but this has almost no effect at solar energies 
where the cross section is essentially due to direct El capture. 

Direct ^Be(p, 7)^B experiments require radioactive targets. It has not been practical to 
use the conventional geometry with large-area, thin targets, and "pencil" beams; instead 
the experimenters were forced to use comparable beam and target sizes. As a result the 
absolute normalization of the cross sections has posed severe experimental problems. 

In the experiments to date, the mean areal density of ^Be atoms seen by the proton beam 
has been determined in one of two ways: 

1. counting the number of ^Be atoms by detecting the 478 keV photons emitted in ^Be 
decay and measuring the target spot size (Wiezorek et al., 1977; Filippone et al, 
1983a,b). 

2. measuring the yield of the ^Li((i,p)^Li reaction on the daughter ^Li atoms that build 
up in the targets as the ^Be decays (Kavanagh, 1960; Parker, 1966, 1968; Kavanagh et 
al, 1969; Vaughn et al, 1970; Filippone et al, 1982). These measurements are made 
on the peak of a broad (P 0.2 MeV) resonance at E^ = 0.78 MeV. 

The first method has the advantage of being direct. The second method has the advantage 
that the ^B produced in the {p, 7) reaction and the ^Li produced in the {d, p) calibration 
reaction can both be detected by counting the beta-delayed alphas, so that detection effi- 
ciency uncertainties largely cancel out. However the second method requires an absolute 
measurement of the total ^Li(d,p)^Li cross section which has turned out to be rather difficult 
to measure correctly. 

The absolute ^Be(p,7)*B cross sections originally quoted from these experiments were 
not consistent with each other, although the shapes of the cross sections as functions of 
bombarding energy were in agreement. Furthermore, the quoted '^L[{d. p)^Li normalization 
cross sections also differed by much more than the quoted uncertainties (values differing by 
up to a factor of two were quoted). However, as pointed out by Barker and Spear (1986), 
even after all the ^Be(p, 7)^B cross sections are renormalized to a common value of the 
^Li((i,p)^Li cross section, the results are not consistent. 

Because poorly understood systematic errors dominated the actTial uncertainties in the 
results, we adopt the following guidelines for evaluating the existing data to arrive at a 
recommended value for 5'i7(0). 

1. We consider only those experiments that were described in sufficient detail that we 
can assess the reliability of the error assignments. 

2. We review experiments that pass the above cuts and make our own assessment of 
the systematic errors, using information given in the original paper plus more recent 
information (such as improved values for the ^Li(o?,p)®Li cross section) when available. 
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The only low-energy ^Be(j9, 7)^B measurement that meets these criteria is the experi- 
ment of Filippone et al. (1983a,b) at Argonne. Filippone et al. (1983a,b) obtained the 
areal density of their target by counting the 478 keV radiation from "^Be decay and also by 
detecting the {d,p) reaction on the ''Li produced in the target by '^Be decay. The Argonne 
experimenters made two independent measurements of the ^Li((i,p)^Li cross section [Elwyn 
et al. (1982) and Filippone et al. (1982)]. These two determinations were consistent. In 
addition, Filippone et al.^s (1982) gamma-ray counting and {d,p) normalization techniques 
gave results in excellent agreement. 



C. The ^Li(d, p)*Li cross section on the E = 0.6 MeV resonance 

Strieder et al. (1996) give a complete listing of existing ^Li{d,p)^Li cross-section mea- 
surements. The results scatter from a maximum value of (211 ± 15) mb (Parker, 1966) to a 
minimum of (110±22) mb (Haight, Matthews, and Bauer, 1985). We obtain a recommended 
value for the ^Li((i, p)®Li cross section by applying the same criteria used above in evaluating 
the ^Be(p, 7)'^B data. The experiments that pass our selection criteria are listed in Table 
p| . The absolute cross sections given in the first three rows of Table |lll| are based on target 



areal densities determined from the energy loss of protons (McClenahan and Segal, 1975) 
or deuterons (Elwyn et al, 1982 and Filippone et al, 1982) in the targets. These results 
therefore share a common systematic uncertainty in the stopping powers. Filippone et al. 
(1982) cite evidence that the tabulated stopping powers were accurate to 5%, but quote an 
overall target thickness uncertainty of 7%. Elywn et al. (1982) quote a ~7.5% uncertainty 
in the stopping power. McClenahan and Segal (1975) quote a target thickness uncertainty 
of 10%. 

The last two entries in Table |T| differ from those given by the authors. The next-to-last 
row was obtained by combining Filippone et a/.'s (1982) two independent, but concordant, 
normalizations of their target thickness. The normalization based on counting the 478 keV 
photon activity from ^Be decay implies a corresponding areal density of ^Li in the target, 
and hence can be used to give an independent absolute normalization to their '^Li{d,p)^Li 
cross section. We obtained the next-to-last value in Table |TTl| by requiring that Filippone et 
a/.'s (1982) measured ''Li + d yield corresponded exactly to their measured '^Li areal density 
inferred by counting the 478 keV photons. Finally, the errors on the '^Li{d,p)^Li cross section 
quoted by Strieder et al. (1996) are unrealistic. Strieder et al. (1996) used a ^Li beam on 
a D2 gas target. They normalized their target density and geometry factor to the ''Li + d 
elastic scattering cross section, which they assumed had reached the Rutherford value at 
their lowest measured energy E = 0.1 MeV. However, their data (see their Fig. 5) do not 
show that the '^Li{d,p) cross section divided by the Rutherford cross section had become 



constant at this energy. Therefore, in the last row in Table |lTl| , we replace Strieder et a/.'s 
(1996) quoted 5% error in the elastic scattering cross section with an 11% uncertainty which 
is the quadratic sum of the 10% uncertainty in the absolute ^Li((i, p)^Li cross section quoted 
by Ford (1964) [Ford's absolute normalization agrees very well with that of Filippone et al. 
(1982)] and a 5% uncertainty in relative normalization of Strieder et a/.'s (1996) data to 
those of Ford. 

We obtain our recommended value for the ''Li((i,p)^Li cross section by the following 
somewhat arbitrary procedure necessitated by the fact that McClenahan and Segal (1975) 
do not give enough information to do otherwise. We assume that each of the first three 
entries in Table had assigned a 7% uncertainty to the stopping power and subtract this 
error in quadrature from the quoted uncertainties. We then combine the resulting values as 
if they were completely independent and then add back a conservative 7% common-mode 
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error. This value is then combined with those of the last two rows in Table |T| which are 
treated as completely independent results. 



D. Indirect experiments 

Two indirect techniques have been proposed that may eventually provide useful quantita- 
tive information on the low-energy ^Be(p, 7)^6 reaction: dissociation of ^B's in the Coulomb 
field of heavy nuclei (Motobayashi et ai, 1994), and measurement of the ^Be + p nu- 

clear vertex constant using single- nucleon transfer reactions (Xu et ai, 1994). Motobayashi 
et al. (1994) quote a "very preliminary value" of S'i7(0) = (16.7 ±3.2) eV b. Measurements 
at low bombarding energies may also provide a constraint of 5*17 (Schwarzenberg et al., 1996; 
Shyam and Thompson, 1997). 

At this point it would be premature to use information from these techniques when 
deriving a recommended value of 5*17(0) because the quantitative validity of the techniques 
has yet to be demonstrated. 

What would constitute a suitable demonstration? In the case of the Coulomb dissociation 
studies, we need a measurement of a dissociation reaction in which radiative capture can also 
be studied directly; the ideal test case will have many features in common with ^Be(p, 7)^B, 
i.e., a low Q-value, a non-resonant El cross section, and similar Coulomb acceleration of 
the reaction products. However, the dissociation cross section has a very different depen- 
dence on the multipolarity than does the radiative capture process. Although ^^0(p,7)^^F, 
^H(a,7)'^Li, and ^^C(p, 7)^^N each has some of the desired properties, a suitable test case 
in which the dominant capture multipolarity is El and the nuclear structure is sufficiently 
simple has not yet been identified. On the other hand, a measurement of the ^'^F — ^^O + p 
vertex constant and the prediction, using the measured vertex constant, of the ^^0(p, 7)^''F 
capture reaction at low energies will provide a good test of the vertex-constant technique. 

To be useful as tests, the indirect calibration reaction and the comparison direct reaction 
must both be measured with an accuracy of 10% or better. Otherwise, one cannot have 
confidence in the method to the accuracy required for the cross section of the ^Be(p, 7)^B 
reaction. 



E. Recommendations and conclusions 

We recommend the value 

5i7(0) = 19^^ eV b , (28) 

where the la error contains our best estimate of the systematic as well as statistical er- 
rors. The recommended value is based entirely on the ^Be(p, 7)*B data of Filippone et al. 
(1983a,b) and is 15% smaller than the previous, widely-used value of 22.4 eV b (Johnson et 
al., 1992) that was based upon a weighted average of all of the available experiments. The 
cross sections were obtained by combining Filippone et al.^s (1982) two independent deter- 
minations of the target areal density [for the ^Li((i, p)^Li method we used the recommended 
cross section in Table |IT1|, and extrapolated these to solar energies using the calculation of 
Johnson et al. (1992). It is important to note that in the region around Ep = 1 MeV where 
the two data sets overlap, Filippone et a/.'s (1983a,b) cross sections agree well with those of 
Vaughn et al. (1970). [We renormalized the Vaughn et al. (1970) data to our recommended 
'^Li{d,p)^Li cross section.] 
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Because history has shown that the uncertainties in determining this cross-section factor 
are dominated by systematic effects, it is difficult to produce a 3cr confidence interval from a 
single acceptable measurement. Instead, we quote a "prudent conservative range," outside 
of which it is unhkely that the "true" 5*17(0) lies 

5i7(0) = 19tl eV b . (29) 

Past experience with measurements of the ^Be(p, 7)^B cross section demonstrates the 
unsatisfactory nature of the existing situation in which the recommended value for S{0) 
depends on a single measurement. It is essential to have additional ^Be(p, 7)^B measure- 
ments, to establish a secure basis for assessing the best estimate and the systematic errors 
for SniO). 

Experiments with ''Be ion beams would be valuable. Such experiments would avoid many 
of the systematic uncertainties that are important in interpreting measurements of proton 
capture on a ^Be target. For example, experiments performed with a radioactive beam 
can measure the beam-target luminosity by observing the recoil protons and Rutherford 
scattering. But the ^Be-beam experiments will have their own set of systematic uncertainties 
that must be understood. Fortunately, experiments with ^Be beams are being initiated at 
several laboratories and results from the first of these measurements may be available within 
a year or two. 

Various theoretical calculations of the ratio of the S*- value at 300 keV and at 20 keV differ 
by several percent. Since these differences will be difficult to measure, yet will contribute 
to the systematic uncertainty in future precise determinations of the solar S-value, a careful 
theoretical study should be made to try to understand the origins of the differences in the 
extrapolations. 



F. Late Breaking News 

In a recent experiment Hammache et al. (1998) measured the cross section at 14 energy 
points between 0.35 and 1.4 MeV (in the center-of-mass system), excluding the 1^ resonance 
energy range. In this experiment two different targets were used with different activities but 
similar results. Hammache et al. determined the ^Be areal density using the two methods 
employed by Filippone et al. (1983a,b) and find consistent results. The measured cross 
section values are in excellent agreement with those of Filippone et al. over the wide energy 
range where both experiments overlap. 

Weissman et al. (1998) report a new measurement of the ''Li((i, p)^Li cross section, 
155±8 mb . The authors also draw attention to the importance of the possible loss of product 
nuclei from the target in cross section measurements performed with high-Z backings. The 
net result of including this new measurement of the ''Li((i,p)*Li cross section together with 
the values given in Table |ITl| , combined with estimates of the effect of loss of product nuclei 
on the previously computed values of 5*17, is a cross-section factor for '^B production that is 
very close to the best-estimate given in Eq. (|29|) . 



IX. NUCLEAR REACTION RATES IN THE CNO CYCLE 

The CNO reactions in the Sun form a polycycle of reactions, among which the main 
CNO-I cycle accounts for 99% of CNO energy production. The contribution of the CNO 
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cycles to the total solar energy output is believed to be small, and, in standard solar models, 
CNO neutrinos account for about 2% of the total neutrino flux. CNO reactions have been 
studied much less extensively than the pp reactions and therefore, in some important cases, 
we are unable to determine reliable error limits for the low-energy cross-section factors. 

Network calculations show that three reactions primarily determine the reaction rates 
of the CNO cycles. The three reactions, "'^^N(p, 7)-'^^0, ^^0(p, 7)"'^'^F, and ^'^0{p,ay'^N, are 
considered in some detail in this review. With a nuclear reaction rate almost 100 times slower 
than the other CNO-I reactions, the reaction ^^N(p, 7)^^0 determines, at solar temperatures, 
the rate of the main CNO cycle. The ^^N and ^^O neutrinos have energies and fluxes 
{El, < 1.8 MeV, (pi^i^CNO) / (j)uCBe) ^ 0.2) comparable to the ^Be neutrinos. The production 
of ^^F neutrinos, with a flux two orders of magnitude smaller, is determined by the reaction 
^^0(j9, 7)^'^F in the second cycle, while ^'^0(p, a)^^N closes the second branch of the CNO 
cycle. 

Figure |^ shows the most important CNO reactions. 



A. i^N(p,7)15q 

1. Current Status and Results 

A number of measurements of the ^^N(p, 7)^^0 cross section have been carried out over 
the past 45 years. Most recently, Schroder et al. (1987), measured the prompt capture 7 
radiations from this reaction at energies as low as Ep = 205 keV; the 1957 measurements of 
the residual /^^-activity of ^^O carried out by Lamb and Hester (1957) between Ep = 100 
and 135 keV remain the lowest proton bombarding energies to be reached in this reaction. 
The solar Gamow peak is at Eq = 26 keV. Three other experiments are available: Hebbard 
and Bailey (1963), Pixley (1957), and Duncan and Perry (1951). 

Table ^ summarizes the measurements and the 5- values determined in previous publi- 
cations, as well as our recommendations. 

As emphasized by Schroder et al. (1987), the relative contributions to the reaction 
mechanism are not fully understood. While Hebbard and Bailey (1963) analyze the data in 
terms of hard-sphere direct-capture mechanisms to the ground, 6.16 MeV, and 6.79 MeV 
states of ^^O, Schroder et al. (1987) find a significant contribution to the ground-state 
capture from the subthreshold resonance at Eji = —504 keV, which corresponds to the 
6.79-MeV state. The agreement of the S'-values recommended by Schroder (1987) and by 
Hebbard and Bailey (1963) seems therefore accidental. The unexplained 40% correction to 
the 7-ray detection efficiency of Schardt, Fowler, and Lauritsen (1952) [an experiment on 
^^N(p, a)^^C used as a cross-section normalization by Hebbard and Bailey (1963)] and the 
anomalous energy dependence of the cross sections in Hebbard and Bailey's (1963) analysis 
argue against inclusion of their results in a modern evaluation of S(0). The lack of a refereed 
publication describing the work of Pixley (1957), and the use of Geiger-counter technology 
in the pioneering experiment of Duncan and Perry (1951), are responsible for our excluding 
these data from the final evaluation. 
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2. Stopping Power Corrections 



The ^^N(p, 7)^^0 cross sections of Lamb and Hester (1957) are important for our un- 
derstanding of the CNO-I cycle, since the data were obtained over an energy range signif- 
icantly closer to the solar Gamow peak (about 30 keV) than other studies of this reaction 
(see Table [TV] ). Lamb and Hester concluded that the S'-factor for this reaction was essen- 
tially constant over the range of proton beam energies from 100 to 135 keV with a value 
S = (2.7 ± 0.2) keV b. Their measurements were carried out using thick TiN targets and 
hence measured yields were integrated over energy as the beam slowed down in the target. 
They assumed a constant stopping power of 2.35 x 10~^° MeV cm^/atom, a good approxi- 
mation at these energies — a recent tabulation (Ziegler, Biersack, and Littmark, 1985) gives 
values of 2.30 x lO'^o MeV cmVatom at 100 keV and 2.22 x 10"^° MeV cmVatom at 135 
keV. In view of the intense proton beams used by Lamb and Hester, there may have been 
significant hydrogen content in their targets, which would increase the molecular stopping 
power by 10% (for TiNH instead of TiN). 



3. Screening Corrections 
Low-energy laboratory fusion cross sections are enhanced by electron screening [see 



Sec. [11 B| and Assenbaum, Langanke, and Rolfs (1987)]. Screening is a significant effect 
at the low energies at which Lamb and Hester (1957) explored the ^^N(p, 7)^^0 reaction. 
Rolfs and Barnes (1990) show that screening effects become negligible for energy ratios 
E/Ue > 1000, where Ue describes the screening potential. This condition is not satisfied 
for the data of Lamb and Hester (1957). Within the adiabatic approximation (Shoppa et 

al, 1993), the screening enhancement can be estimated as f{E) ^ exp |59.6i?~'^/^|, with 

the scattering energy E in keV. (This estimate has only been verified for atomic targets.) 
Screening, and the change in the half-life of ^^O from 120 to 122.2 seconds, are treated as 
corrections, while the considerations related to stopping power are viewed as included in the 
uncertainties quoted by Lamb and Hester. The screening and lifetime corrections reduce by 
8% the S{0) value that otherwise would be inferred from the Lamb and Hester results. 



4. Width of the 6.79 MeV State 

Schroder et al. (1987) made detailed studies of the radiative capture to the bound states 
of ^^O, finding in one case, the ground-state transition, marked evidence for the influence 
of a subthreshold state, the 6.79-MeV level. They were able to observe the capture to this 
state directly, and could thus obtain a proton reduced width. The gamma width, however, 
is not known. Schroder et al. (1987) extracted the gamma width as a fit parameter, finding 
an on-shell width of 6.3 eV. Including the subthreshold state substantially improves the fit 
to the data at energies as high as Ep = 2.5 MeV. However, at the lowest energies for which 
the ground-state transition was measured, the cross section (on the wings of the 278-keV 
resonance) is not well described by the published fitting function. Since gamma-width of the 
6.79 MeV state is not well constrained, the S'-factor for the ground-state transition might in 
principle increase even more rapidly at low energies than found by Schroder et al. (1987), if 
the data at the lowest measured energies were more heavily weighted in the fitting. 

Fortunately, however, there exists a precise measurement of the gamma width of the 
7.30-MeV analog state in ^^N. Moreh, Sellyey, and Vodhanel (1981) find for that state that 
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= 1.08(8) eV, which would imply for the 6.79-MeV state a width of 0.87 eV if analog 
symmetry were perfect. An example is known, however, of a case {A = 13) of an isovector 
El transition that shows considerable departure (more than a factor of two) from analog 
symmetry, but a factor of seven would be surprising. It appears probable, therefore, that 
the width of the 6. 79- Me V state is not significantly larger than that found by Schroder et 
al. (1987). A direct measurement of the gamma width of the 6.79 MeV state would be 
valuable. 



5. Conclusions and Recommended S- Factor for N{p,^y^O 

The experiments of Schroder et al. (1987) and Lamb and Hester (1957) can be used to 
estimate S{0) and its energy derivative. Schroder et al. (1987) provide the only detailed 
data on the reaction mechanism, finding that S rises at lower energies as a result of the 
subthreshold resonance at Eji = —504 keV, while Lamb and Hester (1957) constrain the 
total cross section at the lowest energies. The extent to which the subthreshold resonances 
affect the extrapolation to astrophysical energies is, however, limited by the known width 
of the analog state at 7.30 MeV in ^^N, and, to a degree, by the total cross section from 
Lamb and Hester (1957). The value quoted by Schroder et al. (1987) is therefore likely 
to represent the maximum contribution from a subthreshold state, and cross sections could 
possibly range down to the values found in the absence of the subthreshold resonance. 
There is an uncertainty in the normalization of the two experiments as well, and the overall 
normalization uncertainty is derived as the quadrature of the individual uncertainties. 

The recommended value, 

S{0) = 3.5+°:^ keV b, (30) 

has been obtained by adopting the energy dependences given by Schroder et al. (1987) 
in the presence and the absence of the subthreshold resonance. The energy dependence is 
parameterized in terms of the intercept S'(O) and S"(0) 

S'{0) = -0.008[5(0) - 1.9] b. (31) 

The available data are insufficient to determine S". 

At the mean energy of 120 keV, the data of Lamb and Hester (1957), for which the 
statistical and normalization uncertainty is 12%, have been corrected as described to give 
S'(120) = 2.48 ± 0.31 keV b. For each choice of energy dependence, those data have been 
converted to zero energy and a weighted average formed with the data of Schroder et al. 
(1987), for which the statistical and normalization uncertainty is 17%. The ra-sigma upper 
limits on the average are a quadrature of 3.7 + r;,(0.45) and 3.2 + 72(0.54) keV b; the lower 
limits are a quadrature of 2.5 — ?7,(0.30) and 1.9 — n(0.31) keV b. This prescription, while 
arbitrary, reflects our view that the resonance and no-resonance extrapolations represent a 
total theoretical uncertainty. Hence the recommended "three-sigma" range is 

S{0) = 3.51^:° keV b. (32) 

Figure ^, adapted from Schroder et al. (1987), shows the extant data; the extrapolations 
shown represent the likely range of theoretical uncertainty. Additional uncertainty from 
normalization is not shown in the figure. 

The uncertainty in the ^^N(p, 7)^^0 reaction rate is much larger than previously assumed, 
and produces comparable uncertainties in the calculated CNO neutrino fluxes. On the other 
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hand, the most important calculated solar neutrino fluxes from the p-p cycle are affected 
by at most 1% for a 50% change in the ^^N(p, 7)^^0 reaction rate, as can be seen using the 
logarithmic partial derivatives given by Bahcall (1989). 

New experiments are necessary to improve the understanding of the capture mechanism 
and the cross sections in -^^N(p, 7)-'^^0. 



B. 16o(p,7)17f 

The rate of ^^F neutrino production in the Sun is determined primarily (see Bahcall 
and Ulrich, 1988) by the rate of the ^^0(p, 7)^^F reaction. A number of measurements of 
the ^^0(p, 7)^^F reaction were made between 1949 and the early 70's, and the data are all 
in relatively good agreement. Tanner's (1959) work is consistent with Hester, Pixley and 
Lamb's (1958) lower-energy measurement. Rolfs' (1973) higher-precision work yields the 
value 

5(0) = 9.4± 1.7 keV b . (33) 

No resonance occurs below Ep = 2.5 MeV and a direct capture model describes well the 
data over the entire energy range studied. Since all of the experimental results are consistent 
with each other, Rolfs' (1973) value is adopted. For the latest work on this reaction, see 
Morlock et al. (1997). 



C. ^^0(p,a)i^N 

The ^"^©(p, a)^^N reaction closes the CNO-II branch of the CNO cycles. The S'-factor for 
this reaction has been particularly difficult to measure or predict at solar energies, because 
of the large number of resonances and the difficulty of detecting low-energy alphas. Rolfs 
and Rodney (1975) suggested that a 66 keV resonance may introduce complications arising 
from the interference of the 5604 and 5668 keV energy levels of ^^O. In 1995, an experiment 
at Triangle Universities Nuclear Laboratory (Blackmon et al, 1995) disclosed a resonance 
located between 65 and 75 keV in a comparison of the alpha yields from ^^O and ^^O targets. 
Experiments done by the Bochum group (Berheide et al, 1992), on the other hand, do not 
show evidence for the resonance, and exclude a resonance of the size seen by Blackmon et 
al. (1995), but only on the basis of a smoothly varying background. The proton partial 
width of Blackmon et al. (1995) is Tp = 22tl neV while Berheide et al. (1992) find Fp < 3 
neV. The Bochum group have recently reanalyzed their data, flnding that a different energy 
calibration procedure and choice of background would change their upper limit to 75 neV 
(Trautvetter, 1997). They also have new radiative capture data that indicate an upper 
limit of 38 neV. Landre et al. (1989) measured the proton reduced width in ^''0(^He, (i)^^F, 
but, because the state is weak in proton stripping, uncertainties in the reaction mechanism 
(multi-step and compound nucleus processes) are reflected in the uncertainty; Vp = 71^*^57 
neV. We recommend using the proton width measured by Blackmon et al. (1995), but 
caution the reader that contradictory data have not been revised in the published literature. 

Table summarizes the numerical results. The presence of a near-threshold resonance 
has a signiflcant, but incompletely quantifled, effect on the ^'^0(p, a)^^N cross section at 
solar energies. 
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D. Other CNO Reactions 



We have recomputed the cross-section factors for the ^^C(p, 7)^^N reaction, combin- 
ing the data of Rolfs and Azuma (1974) and Hebbard and Vogl (1960). We find S{0) = 
(1.34 ± 0.21) keV b, S'{0) = 2.6 x lO'^ b, and ^"(0) = 8.3 x 10"^ b/keV. For the reac- 
tion ^^C(p, 7)^^N, we recommend the most recent determination of the S- value reported in 
Table 0, i.e., the values given by King et al. (1994). 

For the ^^N(p, ao)^^C reaction, we have computed the weighted average cross-section 
factor using the results of Redder et al. (1982) and of Zyskind and Parker (1979) [including 
the more accurate measurement by Redder et al. of the cross section at the peak of the 
resonance]. We find a weighted average of 5(0) = (67.5 ± 4) x 10'^ keV b. The cross-section 
derivatives are 5'(0) = 310 b and 5"(0) = 12 b/keV. 

For the reaction ^^0(p, q;)^^N, only an approximate S-value is given since S{E) cannot 
be described by the usual Taylor series and the original analysis by Lorenz-Wirzba et al. 
(1979) determined directly the stellar reaction rates. Wiescher and Kettner (1982) suggest 
a modification of the rate. Very recently, Spyrou et al. (1997) have measured cross sections 
for the ^^F(p, q;)^^0 reaction, but the ^-factor was not determined at energies of interest in 
solar fusion. 



E. Summary of CNO Reactions 



Table |VI| summarizes the most recently published S"- values and derivatives for reactions in 
the solar CNO-cycle. Since the reaction ^^N(p, 7)^^0 is the most important for calculations 
of stellar energy generation and solar neutrino fluxes, it is treated in detail in Table |3 and 
the recommended values for the cross-section factor and its uncertainties are presented in 
Sec. [IX A5|. Other CNO reactions are discussed in Sec. [IX B|, Sec. [IX Q, and Sec. [IX D|. 



F. Recommended New Experiments and Calculations 

Further experimental and theoretical work on the ^^N(p, 7)^^0 reaction is required in 
order to reach the level of accuracy (~ 10%) for the low-energy cross-section factor that is 
needed in stellar evolution calculations. 



1. Low- energy Cross Section 

The cross-section factor for capture directly to the ground state is expected to increase 
steeply at energies below the resonance energy of 278 keV; direct experimental proof of 
this increase is not yet available. Experiments at the Gran Sasso underground laboratory 
(LNGS) using a 1 kg low-level Ge-detector have shown (Balysh et al., 1994) no background 
events in the energy region near = 7.5 MeV over several days of running. A Ge-detector 
arrangement coupled with a 200-kV high-current accelerator at LNGS [LUNA phase II; 
Greife et al, 1994; Fiorentini, Kavanagh, and Rolfs, 1995; LUNA-Collaboration (Arpesella 
et al, 1996)] would allow measurements down to proton energies of 82 keV (corresponding 
to 1 event per day) and could thus confirm or reject the predicted steep increase in S{E) 
for direct captures to the ground state. Still lower energies might be reached by detecting 
the ^^O residual nuclides via their /3^-decay (Ti = 122 s). 
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2. R-matrix Fits and Estimates of the N(p,^)^^ O Cross Section 

Though not fully described, the fit to the ground-state transition in Schroder et al. (1987) 
seems to be based on single Breit-Wigner R-matrix resonances and a direct-capture (DC) 
model added according to a simple prescription not entirely consistent with i?-matrix theory. 
An alternative approach would be to fit the groundstate transition including direct-capture 
and resonant amplitudes following, for example, the description of Barker and Kajino (1991). 
Proper account should be taken of the target thickness. Elastic scattering data of protons 
on ^"^N should be included in the analysis. 



3. Gamma Width Measurement of the 6. 79 Me V State 

Schroder et al. (1987) suggest a large contribution of the sub-threshold state at 6.79 
MeV in ^^O to the ^^N(p, 7)^^0 capture data, and find that the gamma width of that state 
is 6.3 eV. Other experiments yield only an upper limit of 28 fs (F^ > 0.024 eV, Ajzenberg- 
Selove, 1991) for the lifetime of the 6.79-MeV state. Depending upon the actual width, the 
Variant Doppler Shift Attenuation Method (Warburton, Olness, and Lister, 1979; Catford 
et al., 1983), or Coulomb excitation of a ^^O radioactive beam, might yield an independent 
measurement of this width. Data on the Coulomb dissociation of ^^O could also shed light 
on the partial cross sections to the ground state (but not on the total cross section, which 
includes important contributions from capture transitions into ^^O excited states). 



X. DISCUSSION AND CONCLUSIONS 

Table | summarizes our best estimates, and the associated uncertainties, for the low- 
energy cross sections of the most important solar fusion reactions. The considerations that 
led to the tabulated values are discussed in detail in the sections devoted to each reaction. 

Our review of solar fusion reactions has raised a number of questions, some of which we 
have resolved and others of which remain open and must be addressed by future measure- 
ments and calculations. The reader is referred to the specialized sections for a discussion of 
the most important additional research that is required for each of the reactions we discuss. 

Our overall conclusion is that the knowledge of nuclear fusion reactions under solar con- 
ditions is, in general, detailed and accurate and is sufficient for making relatively precise 
predictions of solar neutrino fiuxes from solar model calculations. However, a number of 
important steps still must be taken in order that the full potential of solar neutrino ex- 
periments can be utilized for astronomical purposes and for investigating possible physics 
beyond the minimal standard electroweak model. 

We highlight here four of the most important reactions for which further work is required. 

• The only major reaction that has so far been studied in the region of the Gamow 
energy peak is the ^He(^He, 2p)'^He reaction. A more detailed study of this reaction at low 
energies is required, with special attention to the region between 15 keV and 60 keV. 

• The six measurements of the ^He(a,7)^Be reaction made by direct capture differ by 
about 2.5(7 from the measurements made using activity measurements. Additional precision 
experiments that could clarify the origin of this apparent difference would be very valuable. 
It would also be important to make measurements of the cross section for the ^He(a,7)^Be 
reaction at energies closer to the Gamow peak. 
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• The most important nuclear fusion reaction for interpreting solar neutrino experiments 
is the ''Be(p,7)^B reaction. Unfortunately, among all of the major solar fusion reactions, 
the ^Be(p,7)^B reaction is the least well known experimentally. Additional precise measure- 
ments, particularly at energies below 300 keV, are required in order to understand fully the 
implications of the new set of solar neutrino experiments, Super-Kamiokande, SNO, and 
ICARUS, that will determine the solar '^B neutrino flux with high statistical significance. 

• The ''^^N(p,7)^^0 reaction plays the dominant role in determining the rate of energy 
generation of the CNO cycle, but the rate of this reaction is not well known. The most im- 
portant uncertainties concern the size of the contribution to the total rate of a subthreshold 
state and the absolute normalization of the low-energy cross-section data. New measure- 
ments with modern techniques are required. 
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TABLE I. Best estimate low-energy nuclear reaction cross-section factors and their estimated 
la errors. 



Reaction 


S{0) 


S'{0) 




(keV b) 


(b) 



iH(p,e+z.e)'H 
iH(pe-,z/e)'H 
^He (^He, 2p) "^He 
^He (a, 7) ^Be 
^He (p, e+z^e)^He 
^Be(e-,i/e)'^Li 
^Be {p, 7) ^B 
i%(p,7) 



4.00 (1 ± 0.007+[j:[!?'i') E - 22 
Eq. dl) 
(5.4±0.4)E + 3^ 
0.53 ±0.05 
2.3E - 20 
Eq. dH) 

hO.OO^ 
-0.002 



0^019+0.004 



3.5 



+0.4 



4.48E - 24 



-3.0E - 4 



See Sec. VIII A 



See Sec. IX A 5 



^ Value at the Gamow peak, no derivative required. See text for 5(0), 5' (0). 



TABLE II. Measured values of 534(0). 



534(0) (keV b) 



Reference 



Measurement of capture 7-rays: 

0.47 ±0.05 
0.58 ±0.07 
0.45 ± 0.06 
0.52 ±0.03 
0.47 ±0.04 
0.53 ±0.03 

Measurement of "^Be activity: 

0.535 ± 0.04 
0.63 ± 0.04 
0.56 ±0.03 



Weighted Mean 



Weighted Mean 



Parker and Kavanagh (1963) 

Nagatani, Dwarakanath, and Ashery (1969) 

Krawinkel et al. (1982)^ 

Osborne et al. (1982, 1984) 

Alexander et al. (1984) 

Hilgemeier et al. (1988) 

0.507 ± .016 



Osborne et al. (1982, 1984) 
Robertson et al. (1983) 
Volk et al. (1983) 
0.572 ± .026 



^As extrapolated using the direct-capture model of Tombrello and Parker (1963). 
^As renormalized by Hilgemeier, et al. (1988). 
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TABLE III. '^Li{d,p)^Li cross section (a) at the peak of the 0.6 MeV resonance'' 



Keierence 


a (mbj 


McClenahan and Segal (1975) 


138 ± 20 


Elywn et al. (1982) 


146 ± 13 


Filinnonp pf nl 


1 40 _|_ 10 

A.'-t(J _1_ i.^ 


Fihppone et al. (1982) (Our evaluation, see text) 


146 ± 19 


Strieder et al. (1996) (Our evaluation, see text) 


144 ± 15 


Recommended value 


147 ± 11 



^see also the discussion of Weissman et al. (1998) in Sec. VIII F 



TABLE IV. Cross-section factor, 5'(0), for the reaction ^^N(p, 7)^^0. The proton energies, Ep, 
at which measurements were made are indicated. 



5(0) 
keV b 


MeV 




Reference 


3.20 ±0.54 

3.32 ±0.12 

3.32 

2.75 

3.12 

2.70 


0.2-3.6 

0.2-1.1 
0.100-0.135 




Schroder et al. (1987) 
Bahcah et al. (1982)'^ 
Fowler, Caughlan, and Zimmerman (1975)*^ 
Hebbard and Bailey (1963) 
Caughlan and Fowler (1962)*^ 
Lamb and Hester (1957) 


q c+0.4 

"J'-J-i.e 






Present recommended value 




'^Compilation and evaluation: no 


original experimental data. 




TABLE V. Near threshold resonance widths for -"^"^©(p, a)"'^'*N 


18 F levels (keV) 


5603.4 5604.9 


5673 


Reference 


(eV) 


43 60 


130 


Mak et al. (1980), Silverstein et al. (1961) 


(eV) 


0.5 0.9 


1.4 


Mak et al. (1980), Silverstein et al. (1961) 


Tp (neV) 




71+40 
'^-57 

<3, <75 


Landre at al. (1989) 
Berheide et al. (1992) 



22 



-4 



Blackmon et al. (1995) 
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TABLE VI. Summary of published ^-values and derivatives for CNO reactions. See text for 
details and discussion. When more than one S-value is given, the recommended value is indicated 
in the table. 



Reaction 


Cycle 


5(0) 
keV b 


5'(0) 
b 


S"{0) 
h keV-i 


Reference 


12C(p,7)13n 


I 


1.34 ±0.21 
1.43 

1.24 ±0.15 


2.6E-3 


8.3E-5 


Recommended; this paper 
Rolfs and Azuma (1974) 
Hebbard and Vogl (1960) 




I 


7.6 ± 1 
10.6 ±0.15 

5.7 ±0.8 

S 9 


-7.8E-3 


7.3E-4 


Recommended; King et al. (1994) 
Hester and Lamb (1961) 
Hebbard and Vogl (1960) 

VVUDQUUiy (xliQ PUWici iiyuz) 


IN [p, 7J u 


T 


o t;+0.4 






bet; XdDlc |i V| 




T 


(6.5±0.4)E±4 
(7.5±0.7)E±4 
5.7E±4 


OIU 

351 


1 9 

11 


rLcCOIIlIIlcIlQcQ, tlllo pdptil 

Redder et al. (1982) 
Zyskind and Parker (1979) 
Schardt, Fowler, and Lauritsen (1952) 


i5N(p,ai)i2c 


I 


0.1 






Rolfs (1977) 


15N(p, 7)1^0 


II 


64 ±6 


2.1E-2 


4.1E-3 


Rolfs and Rodney (1974) 


160(p,7)17F 


II 


9.4 ±1.7 


-2.4E-2 


5.7E-5 


Rolfs (1973) 


i70(p,a)i4N 


II 








Brown (1962) (see Table 0) 
Kieser, Azuma, and Jackson (1979) 


i^O(p,7)18f 


III 


12 ±2 






Rolfs (1973) 


i80(p,a)i5N 


III 


~ 4E±4 






Lorenz-Wirzba et al. (1979) 


18o(p,^)19F 


IV 


15.7 ±2.1 


3.4E-4 


-2.4E-6 


Wiescher et al. (1980) 
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FIG. 1. The figure shows the integrand, Upp{r) x Ud{r), of the nuclear matrix element A versus 
radius (fm). The ordinate is given in units of (fm)~^/^. Here Upp{r) and Ud{r) are, respectively, the 
radial wave functions of the p-p initial state and the deuteron final state. The figure (taken from 
Kamionkowski and Bahcall, 1994) displays the integrand calculated assuming five very different 
p-p potentials. In (a) we show the overlap out to a radius of 50 fm, while in (b) we magnify the first 
5 fm. Even drastic changes in the p-p potential result in relatively small changes of the integrand. 

FIG. 2. This figure is adapted from Fig. 9 of Junker et al. (1997), a recent paper by the 
LUNA Collaboration. The measured cross-section factor S{E) for the '^He(^Hc, 2p)^Hc reaction is 
shown and a fit with a screening potential Ue is illustrated. The Gamow peak at the solar central 
temperature is shown in arbitrary units. The data shown here correspond to a bare nucleus value 
at zero energy of 5(0) = 5.4 MeV b and a value at the Gamow peak of 5(Gamow Peak) = 5.3 
MeV b. 

FIG. 3. Comparison of the energy dependence of the direct-capture model calculation (Tombrello 
and Parker, 1963) with the energy dependence of each of the four S34{E) data sets which cover a 
significant energy range. The data sets have been shifted arbitrarily in order to show the comparison 
of the calculation with each data set. 
[Hi88]: (Hilgemeier et al, 1988) 
[Kr82]: (Krawinkel et al., 1982) 
[Os82]: (Osborne et al, 1982) 
[Pa63]: (Parker and Kavanagh, 1963) 

FIG. 4. Model calculations (Tombrello and Parker, 1963) of the fractional contributions of vari- 
ous partial waves and multipolarities to the total (ground state plus first excited state) ^He(a, 7)^Be 
direct-capture cross section factor. 

FIG. 5. CNO reactions summarized in schematic form. The widths of the arrows illustrate the 
significance of the reactions in determining the nuclear fusion rates in the solar CNO cycle. Certain 
"Hot CNO" processes are indicated by dotted lines. 

FIG. 6. Cross sections for ""^^N (p, 7) ""^^O, expressed as S(E), from extant experimental data. The 
data of Lamb and Hester have been corrected as described in the text. The curves represent the 
low-energy extrapolations that would be obtained under the two assumptions of no subthreshold 
resonance (dotted) at Er = —504 keV, and a resonance of the strength considered by Schroder et 
al. (dashed). 
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